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ABSTRACT 

An  analysis  was  done  of  central  Arctic  Ocean  acoustic  data 
to  determine  the  temporal  and  spatial  characteri sti cs  of 
transient  noise  events.  Digital  ambient  noise  data  from 
the  FRAM  IV  experiment  of  April  1982  were  searched  for 
ambient  noise  transients  using  a  detection  program.  The 
time  series  of  the  resulting  detections  were  examined 
visually  to  categorize  each  detection  as  a  transient, 
artifact  or  false  alarm.  The  transient  events  were  located 
in  space  using  time  delays  between  signal  arrival  at 
different  hydrophones.  The  cross  shape  of  the  FRAM  IV 
horizontal  array  permitted  location  in  both  bearing  anc 
range.  The  source  strength  of  each  event  was  calculated 
using  a  simple  dipole  source  model.  Refraction  and 
scatterring  of  the  acoustic  path  in  the  Arctic  Ocean  was 
taken  into  account. 

The  overall  number  of  events  detected,  and  hence  their 
interarrival  times  and  spatial  density,  were  all  affected 
by  the  background  ambient  noise  level.  The  detection 
program  used  the  same  threshold  si gnal — to— noi =e  level  for 
all  data  tapes,  so  when  ambient  noise  levels  were  low  more 
detections  occurred.  The  mean  inter arrival  time  between 
events  was  100  seconds.  The  interarrival  time  fit  a  J 
shaped  gamma  probability  distribution.  The  number  of 
events  detected  per  area  decreased  with  range  from  the 
array  center.  Half  of  the  events  occurred  within  3000 
meters  of  the  array.  In  this  area  there  were  0.3  events 
per  square  kilometer  per  hour.  The  event  population  showed 
no  predominant  angular  dependence.  The  strengths 
calculated  using  the  simple  dipole  model  had  a  mean  o  430 
k H  overal 1  and  260  kN  during  cuiet  times.  Stronger  events 
occurred  during  t  i  mes  with  high  ambient  ncise  levels. 

Theis  Supervisor:  Dr.  Ira  Dyer 

Title:  Professor  of  Ocean  Engineering 
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CHAPTER  1 
INTRODUCTION 

Thesis  Motivation 

This  thesis  investigates  the  spatial  distribution, 
strength  and  rate  of  occurrence  of  law  -frequency  central 
Arctic  noise  events.  During  the  last  30  years  there  has 
been  a  growing  commercial,  military  and  academic  interest 

in  the  Arctic  region. 

In  this  relatively  unexplored  area  there  is  increasing 
evidence  a-f  rich  mineral  and  petroleum  resources.  Research 
into  methods  of  locating  these  assets  and  constructing 
facilities  to  exploit  them  have  received  much  attention. 
These  facilities  must  be  able  to  withstand  the  harsh  Arctic 
surroundings.  The  study  of  Arctic  acoustics  helps  in 
understanding  the  Arctic  environment  and  climate.  It  has 
been  shown  that  there  is  a  direct  correlation  between  10  to 
20  Hz  ambient  sound  pressure  and  environmental  stresses  and 
moments C 1 0 3 .  The  ability  to  use  acoustic  noise  levels  as 
an  environmental  predictor  would  be  a  useful  tool  in  the 
protection  of  commercial  Arctic  facilities. 

The  Arctic  ocean  serves  as  a  military  arena  for 
several  submarine  fleets.  The  underice  environment  makes 
detection  difficult,  increasing  the  stategic  role  of  those 
fleets.  Because  of  the  sound  velocity  profile  of  the 
central  Arctic  there  is  a  surface  duct  which  channels  sound 
for  long  distances.  But,  the  underice  profile  scatters 
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sound  energy,  effectively  filtering  out  high 


f requenciesC 1 1 3 .  The  result  is  that  only  low  frequency 
signals  travel  far  in  the  Arctic,  and  therefore  the  low 
frequency  range  is  the  best  for  detecting  adversary 
submarines.  The  importance  of  understanding  the  low 
frequency  ambient  noise  field  becomes  apparent.  The  actual 
central  Arctic  ambient  noise  level  is  at  times  much  quieter 
than  the  open  ocean,  but  it  contains  unpredictable 
transient  noise  events  which  interfere  with  conventional 
detection  schemes.  It  has  been  hypothesized  that  the 
background  ambient  noise  is  the  summation  of  these 
transients  from  throughout  the  Arctic  basinC73.  Analysis 
of  the  spatial  and  temporal  distribution  of-  these 
transients  is  a  logical  next  step  in  understanding  low 
frequency  noise,  and  improving  our  submarine  detection 
capabi 1 i ty . 

The  academic  challenge  of  the  Arctic  lies  in  the 
sparseness  of  field  data.  The  Arctic  cannot  be  casually 
sampled.  Even  simple  experiments  require  expensive 
expeditions.  The  harsh  environment  takes  its  toll  on 
researchers  and  equipment,  and  reduces  the  amount  of  usable 
data.  Hence,  the  study  of  the  Arctic  is  like  a  jigsaw 
puzzle  with  few  pieces  present.  The  total  picture  remains 
a  stimulating  mystery. 

This  study  analyzes  data  collected  during  the  FRAM  IV 
experiment  by  Massachusetts  Institute  of  Technology  and 


Moods  Hole  Oceanographic  Institute  personnel. 


The  FRAM  IV 


ice  camp  was  located  in  the  Barents  (Nansen)  Abyssal  Plain 
at  approximately  84°  N  by  15°  E,  as  shown  in  Figure  1—1. 
The  ice  was  3  meter  thick  multi-year  pack  ice.  The  ice 
activity  was  low;  there  was  no  ice  ridging  or  lead 
formation  around  the  camp  during  the  experiment. 

The  FRAM  IV  ice  camp  was  set  up  -from  25  March  to  11 
May  1982.  This  study  analyzes  data  taken  between  March 
27th  and  April  22nd.  The  weather  was  mild,  with 
temperature  ranged  from  —35°  to  —4°  C,  and  wind  speed  from 
1  to  23  knots. 

The  ambient  noise  was  sensed  with  a  large  horizontal 
hydrophone  array  which  consisted  of  two  non-uni forml y 
spaced  line  arrays,  crossing  at  right  angles.  The  data 
were  digitally  recorded  on  a  multichannel  system. 
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Figure  1-1  Location  of  the  FRAM  IV  Arctic  exoenment 
conducted  in  the  spring  o-f  1982.  C  1  1  2 
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Thesis  Contents 


My  work  in  this  area  began  long  a-fter  the  FRAM  IV 
ambient  noise  data  were  collected.  The  first  step  was 
finding  the  events  in  the  raw  data.  The  data  consisted  of 
magnetic  tapes  each  containing  20  minutes  of  digitized 
noise  levels.  A  program  was  written  which  searched  the 
ambient  noise  tapes  for  passible  events.  Chapter  2 
di scribes  the  automated  and  manual  techniques  used  to 
accomplish  this  detection. 

These  events  were  then  located  in  space  using  the 
difference  in  arrival  time  between  hydrophones.  This  was 
also  done  with  a  computer  program.  The  program  plotted  the 
arrival  time  delays  against  range  to  a  trial  location,  did 
a  least  squares  fit,  and  chose  the  location  with  the  best 
fit.  This  is  covered  in  Chapter  3. 

The  peak  voltages  for  each  event  were  used  with 
the  dipole  source  model  to  predict  peak  source  strength. 

The  background  ambient  noise  strength  was  also  determined. 
These  strength  calculations  are  found  in  Chapter  4. 

In  Chapter  5  the  distribution  of  event  interarrival 
time  was  determined.  The  event  locations  and  strengths 
were  analyzed,  and  a  spatial  density  found. 

Chapter  6  summarizes  the  key  results  of  this  study. 
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CHAPTER  2 

DETECTION  OF  NOISE  EVENTS 

Twenty  nine  FRAM  IV  ambient  noise  tapes  were  searched 
in  order  to  find  a  population  of  noise  events  for  this 
study.  The  specific  tapes  were  chosen  from  the  possible  67 
in  order  to  cover  the  entire  range  of  days  of  the  FRAM  IV 
experiment.  However,  there  were  several  davs  when  no 
ambient  noise  tapes  were  recorded.  To  help  fill  these  gaps 
five  reverberat i on  tapes  were  also  searched.  These  tapes 
were  recorded  prior  to  the  reverberat i on  shot  being  fired, 
or  they  were  recorded  so  late  in  the  experiment  (SO  minutes 
after  the  shot)  that  reverberati ons  were  no  longer  present. 

The  FRAM  IV  experiemnt  used  the  horizontal  array  of 
omnidirectional  hydrophones  pictured  in  Figure  2-1.  The 
hydrophones  were  suspended  from  the  ice  into  the  water  to  a 
depth  of  93  meters  below  the  air/ice  interface.  The  two 
crossed  lines  of  the  array  allowed  the  possibility  of 
localizing  events  in  space.  Although  26  hydrophones  are 
shown  in  Figure  2-1,  only  24  at  a  time  could  be  used  to 
record  data.  In  most  cases  a  few  of  the  recording  channels 
were  used  for  other  sensors  (geophones  or  hydrophones  used 
in  a  vertical  arrav).  Most  of  the  time  19  to  21  horizontal 
arra.-  h.drr  phone  data  -jere  recorded. 

The  FRAM  IV  ambient  noise  tapes  were  recorded 
digitally.  Figure  2-2  shows  a  schematic  of  the  svstem  used 


FRAM  IV  DATA  ACQUISITION 


HORIZONTAL 

ARRAY 


Figure  2-2  Schematic  of  the  recording  system  used  for  FRAM 
IV  data  collection.  C43 


■for  collecting  the  data.  The  input  from  each  hydrophone 
was  put  through  a  gain  ranging  ampli-fier  and  a  low  pass  80 
Hz  filter.  It  was  sampled  at  250  Hz  and  recorded  on  20 
minute  magnetic  tapeC41.  The  24  channel  recorder  had  a  120 
dB  dynamic  rangeC12]. 


Event  Detector  Program 

The  event  detector  prograta  was  written  to  tare  digital 

data  from  a  FRAM  IV  tape  and  determine  where  in  that  tape 
noise  events  occurred.  The  program  was  originally  written 
to  take  data  directly  from  a  tape  drive,  but  subsequently 
modified  to  take  the  data  from  a  file.  The  framread 
program,  with  a  -head  switch  is  used  to  read  the  tape  into 
the  file.  This  will  eliminate  any  headers  and  then  read 
the  digital  data  straight  into  a  file.  A  FRAM  II  tape  mav 
be  read  into  a  similar  file  using  framread  and  the  switches 
-head  and  — fram2.  The  framread  program  was  written  by 
G.  Duckworth,  and  is  available  to  the  Arctic  Acoustics 
Program  at  MIT. 

The  event  detection  programs  source  codes,  flow  chart, 
and  a  short  users  guide  are  found  in  Appendix  A.  The  event 
detection  program  which  reads  from  a  file  is  called 
hdetect .  The  detection  program  reads  from  a  tape  drive. 

Both  programs  are  written  in  the  c  programming  language  for 
a  UNIX  operating  system. 

The  event  detection  program  follows  the  block  diagram 


/VVV>'^ 


of  Figure  2—3. 


The  program  initialization  portion  defines 


variables  and  constants,  zeros  flags,  and  requests  user 
input  such  as  tape  number,  date,  time  and 

channel s (hydrophones)  to  be  used,  as  well  as,  the  name  of 
input  and  output  files.  After  this  information  is 
requested  from  the  user,  the  program  no  longer  requires 
attent i on . 

The  event  detection  program  then  reads  in  a  file  of 

data,  filters  the  data,  squares  each  data  point,  and  takes 
the  square  root.  The  filter  was  a  Par ks-McCl el  1  an  digital 
20  to  80  Hz  bandpass  filter.  Its  frequency  response  is 
shown  in  Figure  2-4.  The  range  of  this  filter  was  chosen 
to  avoid  the  NyquisJ  Crequency  (125  Hz)  and  hydrophone 
cable  strum  (1-20  Hz),  and  to  be  compatible  with  the  analog 
80  Hz  low  pass  filter  the  data  went  through  before  being 
recorded.  The  data  were  squared  and  then  rooted  to  ensure 
positive  peak  values  far  all  data  points. 

The  lent  portion  of  the  program  used  a  threshold 
detection  scheme  to  check  each  channel  for  possible  noise 
events.  For  a  particular  channel  a  short  average  of  the 
four  most  recent  data  points  was  compared  to  a  long  average 
of  64  recent  data  points.  If  the  ratio  of  short  average  to 
long  average  was  over  a  certain  value  that  channel  would  be 
flagged  for  a  possible  event.  The  time  of  the  flag  and  the 
value  of  the  short  average  were  also  recorded.  All  other 


channels  averages  were  taken  similarly. 


If  at  least  50?.  of 


HI  4j 


*2-3  Block  diagram  a-f  major  modulues  of  the 
ion  program. 
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the  channels  were  -flagged,  an  event  would  be  declared. 

Then  the  program  would  shift  to  the  next  time  increment  of 
data  and  the  process  would  be  repeated. 

A  more  detailed  diagram  of  the  event  detection  module 
is  seen  in  Figure  2-5.  There  are  four  submodules:  reset 
flag,  set  flag,  new  event,  and  deactivate  old  event 
modules.  The  reset  and  set  flag  modules  deal  with  the 
channel  flags  which  trip  when  a  particular  channel 
experiences  a  large  si gnal -to-noi se  ratio  (i.e.  the  RATIO 
of  short  average  to  long  average  exceeds  a  certain  level). 
The  new  event  and  deactivate  old  event  modules  deal  with  an 
active  event  matrix  which  identifies  active  events,  and 
stares  channel  flag  time  and  amplitude  for  each  declared 
event . 

The  reset  flag  module  resets  the  channel  flag  if  it 
has  been  more  than  0.3  seconds  since  the  channel  tripped. 
Spurious  peaks  on  a  channel  might  flag  a  channel 
prematurely.  This  reset  module  prevents  a  number  o-f 
channels  with  spurious  peaks  over  a  long  time  period  from 
being  falsely  declared  an  event.  The  value  of  this  RESET 
DELAY  was  determined  by  examining  known  events  and  noting 
that  about  half  the  channels  tripped  within  a  0.3  second 
peri  ad . 

The  set  flag  module  determines  if  the  short  averane  to 
long  average  RATIO  has  been  exceeded  and,  if  it  has,  the 
module  t)  checks  to  see  if  the  channel  flag  is  already 


tripped.  2)  checks  to  see  1  f  this  detection  is  part  at 
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recently  declared  event,  and  3)  adds  its  detection 
information  to  the  active  event  or  the  channel  flag,  or 
discards  the  information  depending  on  the  circumstances. 

The  short  average  length  o-f  4  <=  0.016  sec)  was  chosen  sc 
that  it  matched  the  length  o-f  the  signal  (0.02  sec).  This 
provides  the  maximum  signal  level  since  this  is  long  enough 
to  get  ail  of  the  signal  and  short  enough  not  so  a /a  rage  :  h 
with  the  lower  surrounding  background  noise  level.  The 
long  average  length  o-f  64  was  chosen  because  a  length  ratio 
o-f  15  to  1  had  been  suggested  by  Ke 1 1  v C 9 3  -for  the  Large 
Aperture  Seizmic  Array  (LASA) .  The  detection  and 
localization  schemes  used  by  this  large  horizontal  arrav 
were  directly  appiicaole  to  the  FRAM  IV  hydrophone  arrav 
data.  The  choice  o-f  detection  RATIO  was  done  through  a 
series  o-f  tests,  and  the  selection  was  made  bv  balancinc 
detection  rate  and  false  alarm  rate. 

when  an  event  is  declared  the  information  in  tne 
channel  flags  is  transferred  to  the  active  event  matrix, 
and  the  channel  flags  are  cleared.  So  if  a  channel  has  a 
detection  and  its  channel  flag  is  not  already  tripped,  the 
set  flag  module  must  first  see  if  the  detection  belongs  to 
a  recently  declared  event.  If  the  active  event  already  has 
1  r  i  t  c t.  r.  f  acgea  .  the  inf  o  rmat  l  on  is  re  n  j  *ceo  c  v  f  h  a 
new  detection  on l v  if  their  times  differ  bv  less  than  0.  i2 
seconds  and  the  new  detection  amDlitude  is  cr~eater.  This 


EVENT  DELAY  time  of  0.5  second  was  chosen  because 
inspection  o-f  the  known  events  revealed  that  most  channe'  5 
tripped  within  a  0.8  second  period.  The  RESET  DELAY  (0.03 
sec)  plus  the  EVENT  DELAY  result  in  this  0.8  second  look 
-for  each  event. 

The  output  o-f  the  event  detection  program  is  an  output 
-file  which  contains  the  time  each  event  was  declared,  which 
channels  were  flagged  and  the  time  del av  and  peal  amplitude 
for  each  channel.  The  time  delays  are  relative  to  the 
earliest  channel  flag  time,  so  one  of  the  channels  always 
has  a  zero  time  delay. 

The  first  version  of  the  detection  program  was  written 
to  take  a  new  short  average  and  long  average  at  everv  data 
point  (every  0.004  seconds).  This  program  took  4  to  3 
hours  to  search  a  20  minute  data  tape.  A  concession  to 
speed  was  made  and  the  program  changed  to  compute  averages 
at  every  fourth  sample  point  (every  0.016  seconds).  This 
reduced  the  accuracy  of  the  time  delays  and  the  abilitv  of 
the  program  to  pick  up  events.  The  RATIO  had  to  be  lowered 
in  order  to  get  the  same  detections  which  were  obtained 
previ ousl y . 

Studies  to  find  the  best  si gnal -to-noi se  RATIO  were 
conducted  several  times.  Development  of  the  LASA  detection 

system  had  revealed  that  a  7  dB  si  gnal -to-noi  se  rati  d  s 

needed  for  75"/.  detect  1  on  C  8  3 .  This  equates  to  a  5  to  1 
ratio  of  signal  power  to  noise  power.  Since  I  was  working 


s' 


with  pressure  vice  power  I  used  2.2  as  my  starting  RATIO. 


This  RATIO  engul-fed  me  in  -false  alarms.  A  quick  study  was 
done  around  the  2.4  level.  The  first  10  minutes  of  tape 
4013  were  run  at  RATIOS  of  2.4,  2.45  and  2.48.  This  tape 
had  been  visually  examined  in  detail  previously,  so  the 
events  were  known.  The  results  are  shown  in  Table  2-1. 
Also  shown  in  Table  2-1  are  the  results  of  a  second  study, 
done  after  the  program  had  been  changed  to  average  lees 
often . 


Table  2-1  Determination  of  the  Best  RATIO 


Average  taken 

at  every  data  point 

RATIO 

detection 

2.4 

947. 

2.45 

757. 

2.48 

7 1  7. 

false  alarm  rate 
257. 

137. 

87. 


Average  taker,  at  every  fourth  data  point 
RATIO  detection  rate 

2.3  767. 

2.38  717. 

2.4  597. 

2.5  597. 


false  alarm  rate 
467. 

297. 


oetec  l n  r  ct e  —  4  event  dt-t&ct  i^n^  4  or  inown  e\  n  t  =? 

FA  rate  =  4  non-event  detections  /  total  detections 


Notice  how  the  detection  rate  has  decreased  and  the 


false  alarm  rate  increased  as  a  result  of  only  averaging  at 
every  fourth  data  point.  Averaging  less  frequently  means 
there  is  a  smaller  probability  that  the  data  points  to  be 
averaged  will  all  lie  near  the  peak  amplitude  of  the 
signal.  The  signal  level  is  generally  lower  than  that 
detected  when  averaging  every  data  point,  and  a  lower 
s i gna 1  -t a— no i se  RATIO  must  be  usee  to  detect  the  same 
events.  But  when  the  si gnal -to— noi se  RATIO  is  lowered  the 
false  alarm  rate  increases. 

The  RATIO  of  2.38  was  settled  an.  This  is  a 
compromise  which  gives  a  detection  rate  which  finds  most 
high  and  medium  strength  events,  and  which  has  a  tolerable 
false  alarm  rate.  Because  the  detection  rate  is  less  than 
1 QQV.  (71X)  there  were  events  present  which  could  be  seen 
visually,  but  were  not  aided  up  bv  the  detection  orogr-m. 
The  RATIO  could  have  been  adjusted  to  detect  all  events 
seen  visually,  but  at  the  cost  of  a  multitude  of  false 
alarms.  The  RATIO  was  kept  at  2.38  and  used  for  the 
detection  of  all  data  tapes. 

The  final  version  of  the  hdetect  program  read  digital 
data  from  a  framread  file,  detected  passible  events  using 
the  less  frequent  averaging  scheme,  and  supplied  the  even- 
riie  nc  channel  ti  ne  del  a  vs  and  amoi  i  tunes  t  -  an  o  -  d  •- 
file.  Once  the  RATIO  had  been  sat i s f actor  1 1 y  set  the 
program  was.  used  to  search  the  FRAM  IV  ambient  noise  tares 


.  r.  /  v. 


V  -*  V  /  '/  V/  V 


.*  >  •  "j*  V 


Lyftjij '  y*ij 


>s  are  usually  within  3  seconds  of  each  ether  . 
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been  off  by  as  much  as  13  seconds  in  one  case.  The  time 
difference  between  the  two  methods  is  greater  at  the  end  of 
a  tape,  and  is  likely  to  occur  after  a  particularly  strong 
event  has  taken  place  (though  there  were  times  when  time 
di screpancies  developed  without  strong  events  present.  and 
also  many  strong  events  existed  which  did  not  induce 
discrepancies).  Typically,  there  might  be  no  time 
difference  at  the  *.rst  part  cf  the  tape,  then  after  a 
strong  event  a  three  second  discrepancy  would  De  seen  ant 
tliis  would  be  consistent  until  the  end  of  the  tape, 
because  the  errors  did  not  appear  randomly  throughout  the 
tape,  and  because  they  developed  impulsively.  I  believe 
that  the  problem  lies  in  the  time  counter  of  the  event 
detection  program  becoming  offset  from  the  time  of  the  raw 
data,  perhaps  because  of  short  r?ccrdc  in  the  raw  da*- a. 

time  keeper  in  the  e  .  en  t  of  a  short  r-e-mr  d .  f.  dr-act  l  nc, 
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for  events,  and  no  further  program  development  was  done. 
Weaknesses  in  the  program  were  subsequently  discovered,  but 
have  not  been  corrected. 

The  biggest  problem  is  the  accuracy  of  the  event  time 
(the  time  when  50’/.  of  the  channels  have  been  flagged).  An 
event  time  reported  by  the  event  detection  program  will  not 
exactly  match  that  found  by  plotting  the  time  series.  The 


f  hi  s  fiflv  eliminate  t  l  me  d  i  a.-  r  “<->  ^nr  ,  or  <--p  ]  e>rn  . 


It  has  already  been  mentioned  that  the  accuracy  of  the 
time  delays  deteriorated  when  the  program  was  changed  to 
run  mare  quickly.  This  also  made  the  detection  o-f  the 
event  peak  amplitude  less  likely.  As  a  result,  in  order  to 
get  accurate  locations  and  source  strengths,  both  time 
delays  and  peak  voltages  had  to  be  taken  manually  from  time 
series  plots  o-f  each  event. 

The  other  major  problem  o-f  the  event  detector  program 
is  that  it  does  not  discriminate  between  an  Arctic  noise 
transient  and  an  artifact  such  as  an  air  gun  blast  or  a 
reverberati on  shot.  Short,  strong  signals  are  reported  as 
passible  events.  Adding  this  discrimination  to  the  program 
is  the  next  step  in  improving  its  usefulness. 

Visual  Confirmation 

Visual  confirmation  was  required  for  all  possible 
noise  events  in  order  to  eliminate  artifacts,  false  alarms 
and  multiple  detections  of  the  same  event.  In  addition,  in 
a  few  cases  visual  confirmation  revealed  two  events  where 
there  had  only  been  one  detection. 

The  event  detection  program  was  designed  to  preclude 
the  need  for  platting  a  time  series  of  each  event.  The 
output  of  the  program  contains  time  delay  and  voltage 
amplitude  information  which  can  be  used  directly  in  the 
location  program.  However,  because  of  the  decreased 
accuracy  of  the  time  delay  and  amplitude  information,  and 


because  of  the  event  time  discrepancy  mentioned  previously, 
it  was  necessary  to  plot  the  time  series  o-f  each  event. 

The  -first  step  o-f  the  visual  confirmation  is  to  review 
the  tape  log  for  any  artifacts  that  may  have  occurred 
during  the  recording.  The  times  are  noted,  and  these  are 
compared  to  the  event  times  given  by  the  detection  program. 
Then  a  time  series  of  the  artifact  was  plotted  to  determine 
which  detections  were  associated  with  it.  In  general,  an 
artifact  such  as  an  air  gun  blast  did  not  affect  detections 
for  over  20  seconds. 

The  visual  confirmation  portion  of  the  procedure 
evolved  from  a  very  limited  look  only  at  events  which  could 
not  be  located  with  the  detection  program  generated  time 
delays,  to  a  three  step  plotting  procedure  for  each  event. 
During  the  early  period  of  this  work  the  hdetect  program 
output  was  used  directly  as  the  input  to  the  location 
program.  The  location  program  used  the  time  delays  to 
determine  the  event's  location  in  space.  Those  events 
which  could  not  be  located  needed  a  closer  look,  and  so 
their  time  series  were  plotted.  The  plots  were  made  of  a  2 
second  period  including  the  event  time  given  by  the  hdetect 
program.  Often  there  was  no  apparent  event  in  this  time 
series  plot,  and  the  detection  was  declared  a  false  alarm. 
When  an  event  did  plot,  manual  time  delays  were  taken  and 
used  to  locate  the  event.  These  manual  delays  located 
these  events  with  better  accuracy  then  the  hdetect 


generated  time  delays. 
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It  soon  became  apparent  that  the 


best  answers  would  be  obtained  by  taking  manual  time  delays 
o-f  al  1  events.  Trying  to  localize  the  events  with  the 
program  generated  time  delays  was  dropped  -from  the 
procedure,  and  the  -first  step  after  getting  the  hdetect 
program  output  became  doing  a  2  second  time  series  plot  o-f 


each  event. 


-f  ter 


iozen  tapes  had  been  analyzed  in  this  name 


the  discovery  o-f  the  event  time  discrepancy  was  made. 

A-fter  platting  a  dozen  events  right  at  the  time  shown  by 
the  program,  the  final  two  dozen  event  o-f  tape  4009  all 
appeared  to  be  -false  alarms.  The  quality  of  the  tape  was 
good  (low  background  noise),  so  this  seemed  highly 
suspicious.  A  broader  search  of  the  time  around  each  event 
showed  that  the  final  two  dozen  events  were  not  false 
alarms,  but  were  events  with  times  5  seconds  different  than 
those  indicated  by  the  program,  so  that  none  of  those 
events  had  shown  up  on  the  2  second  time  series  plats.  The 
method  of  visually  confirming  events  was  changed  so  that  a 
waterfall  plot  was  made  around  the  time  of  each  event. 

This  showed  the  exact  time  of  the  event,  and  helped  discern 
the  pattern  of  time  discrepancy  between  time  series  and  the 
hdetect  program.  Once  the  pattern  was  found  events  were 
easy  to  * i nd  and  false  al arms  could  be  noted.  The 
confirmed  events  were  then  plotted  with  2  second  time 
ser i es . 


<^.1 


i  he  -final  step  at  visual  confirmation  was  simply 
platting  and  replatting  the  events  to  the  proper  nair  sc 
that  the  higher  voltage  amplitudes  would  not  be  cut  off. 
These  peak  voltages  were  manually  taken  from  the  time 
series  plots  and  used  to  determine  source  strength. 

The  final  method  used  for  visual  confirmation  was: 

1)  Check  tape  log  for  artifacts,  and  eliminate 
those  from  further  analysis. 

2)  Plot  a  waterfall  time  series  around  each 
possible  event,  separate  real  events  from  false  alarms,  and 
find  true  event  times. 

3)  Plot  2  second  time  series  of  each  real  event, 
adjusting  gain  to  keep  from  clipping  higher  voltages. 

Using  this  technique  certainly  reduced  the  false  alarm 
rate.  A  breakdown  of  the  detection  statistics  for  tapes 
that  had  been  examined  by  both  methods  is  found  in 
Table  2-2. 


Table  2-2  Detection  Statistics  for  Two  Visual 
Confirmation  Methods 


Art i f  acts 


Events 


False  Alarms 


□rignal  Method 
w/  2  sec  plots 


1 6 . 97. 


37 . 57. 


45 . 67. 


Ultimate  Method 
w/  waterfall  plot 
2  sec  c 1  at 


16.27. 


65. 47. 


IB.  4V. 


( h'er  cen  t  ages  of  detections  classed  in  each  category) 


The  human  interpreter  was  a  necessary  tool  in  this 
scheme.  There  was  not  necessarily  a  one-to-one 
correspondence  between  events  and  detections.  There  were 
cases  where  a  strong  event  would  cause  multiple  detections, 
and  cases  where  two  events  occurred  at  the  same  time  and 
caused  only  a  single  detection.  In  some  cases  a  series  of 
detections  seemed  to  be  an  event  and  an  echo,  or  perhaps 
straf  .  This  would  be  counted  as  a  single  event. 

The  method  o-f  determining  whether  a  detection  was  a 
-false  alarm  or  a  weak  event  was  sometimes  difficult.  In 
general,  if  the  detection  program  indicated  a  possible 
event,  "something"  could  be  seen  on  the  waterfall  plot. 

The  detection  was  dismissed  as  a  false  alarm  if  no  pattern 
for  taking  time  delays  could  be  seen.  (Because  of  the 
shape  of  the  hydrophone  array  there  were  consistent 
patterns  of  time  delays  depending  upon  the  direction  to  the 
event.)  Presumably  the  false  alarm  rate  depends  upon  the 
training  and  attention  of  the  human  interpreter. 

Manual  time  delays  were  taken  from  an  arbitrary 
reference  to  the  crossing  of  the  largest  peak  to  peak 
amplitudes,  as  shown  in  Figure  2-6.  For  most  events  this 
was  clear,  but  for  weak  or  complex  events  some  intuition 
was  needed. 

Voltage  amplitudes  were  tal en  as  the  maximum  peak 
voltage  in  the  event  signal.  All  were  taken  as  magnitudes 


Figure  2-6  Sketch  showing  point  at  measurement  tor  event 
time  delays.  The  event  is  timed  at  its  zero  crossinq 
between  the  largest  pair  at  positive  and  negative  peaks. 
(The  measurement/anal ysi s  system  has  a  polarity  ot  negativ* 
voltage  tor  positive  pressure. ) 
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All  types  of  noise  event  signatures  previously 
observed  by  DyerC7j  were  seen  in  the  ambient  noise  tapes  I 
evaluated.  The  majority  of  events  were  pops  and  extended 
pops.  There  were  also  a  few  whines  and  straf  events. 

While  signature  types  were  noted  in  general,  the  signature 
type  of  each  individual  event  was  not  recorded. 
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I  have  assumed  that  the  strangest  peak  pressure  sensed 
at  the  hydrophone  is  due  to  a  waterborne  acoustic 
propagation  path  -from  a  source  located  in  the  ice  sheet. 

The  signal  enters  the  water  near  the  source  and  propagates 
directly  toward  the  hydrophone. 

I  have  assumed  that  the  signal  does  not  bounce  o-f-f  the 
ocean  bottom  or  the  ice  canopy  before  reaching  the 
hydrophone.  Paths  bouncing  off  the  bottom  would  produce 
signals  with  much  lower  energy  than  the  direct  path  signal, 
and  can  be  ignored.  Signals  bouncing  off  the  ice  canopy 
are  too  energetic  to  ignore  but,  as  I  show  subsequently, 
they  do  not  affect  the  time  delay  computations  significantly. 

The  location  program  is  based  on  arrival  times  being 
related  to  slant  range,  R  ,  and  does  not  take  the  upward 
refraction  of  the  acoustic  path  into  account.  The  impact 
this  has  on  the  results  is  discussed  in  the  next  section. 

The  location  program  takes  as  input  a  file  of  time 
delays  and  voltage  amplitudes,  and  outputs  a  file 
containing  the  best  event  location,  sound  speed,  and 
standard  deviation.  It  also  computes  source  strength  based 
on  the  voltage  amplitude  inputs,  the  event  location  and  a 
spherical  spreading  loss.  This  feature  was  originally 
included  so  that  the  source  strength  could  be  computed 
directly  from  the  event  detection  program  outputs.  Since 
the  peak  voltages  recorded  by  the  detection  program  are  not 
as  accurate  as  those  done  by  hand,  and  since  the 
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transmission  loss  does  not  -follow  simple  spherical 


spreading,  these  computed  source  strengths  were  not  used 
for  any  part  of  this  study. 

The  FQUAK  program  set  up  a  grid  of  points  around  a 
specified  center  position.  The  grid  consisted  of  a  point 
every  lOO  meters  from  -5000  to  +5000  meters  in  both  the 
and  y  directions.  This  resulted  in  100  >:  100  test 

locations.  When  the  best  test  location  was  found  the 
interval  spacing  was  reduced  to  every  10  meters,  and 
another  10,000  test  locations  were  generated  using  the  best 
location  of  the  first  round  as  the  new  center.  The  process 
was  repeated  with  a  1  meter  interval  to  get  the  final 
answer.  The  scheme  evaluated  a  total  of  30,000  test 
locations,  covered  a  range  out  to  5000  meters,  and  cook 
about  20  minutes  to  run. 

I  noted  that  a  significant  number  of  events  found  with 
FAiJAK  were  at  the  range  limit  of  5000  meters.  The  program 
location  was  written  to  search  a  larger  area  faster.  The 
fineness  of  the  grid  was  decreased  to  20  20  vice  100  :: 

100.  A  1000  meter  interval  was  added  to  enable  the  program 
to  search  out  to  10,000  meters.  This  reduced  the  total 
number  of  test  locations  to  1600  (20  >:  20  >: 4 )  ,  and  the  time 
to  one  minute.  location  gave  answers  which  were  very 
consistent  with  F0UAK,  except  in  one  particular  s;  teat  l  o-  . 

The  wider  grid  size  led  to  one  problem.  The  location 
program  sometimes  found  the  lowest  standard  deviation  for  a 
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point  in  the  quadrant  directly  opposite  the  true  location. 
This  was  suggested  by  the  sound  speed  being  reoorted  as 
approximately  -1440  m/sec,  as  illustrated  in  Figure  3—3. 
This  problem  was  solved  by  modifying  location  to  make  the 
program  f  ine locate .  This  program  used  the  grid  size  and 
spacing  of  FQUAK ,  and  centered  the  search  so  that  the  user 
could  designate  which  of  the  quadrants  would  be  searched. 

A  casual  look  at  the  manual  time  delays  of  an  event  eas: 1 v 
reveals  the  appropriate  quadrant.  This  program  works  well, 
but  is  as  slow  as  the  original  FQUAK.  It  was  used  rarely. 

As  with  the  original  FQUAK ,  I  began  to  notice  that 
some  events  were  located  at  the  range  limit  of  the  location 
program.  This  led  to  the  modification  of  the  location 
program  to  form  the  program  far  1 ocate .  This  program  uses 
the  location  grid  size  and  fineness,  but  allows  the  center 
point  to  be  any  of  the  far  corners  of  ihe  original  location 
grid,  or  at  the  limit  range  at  each  of  the  cardinal  points. 
This  is  shown  in  Figure  3—4.  This  allowed  events  to  be 
located  out  to  20,000  meters. 

The  location  orogram  source  code  and  a  brief  user's 
manual  are  found  in  Appendix  B.  This  program  was  written 
in  the  c  programming  language  for  the  UNIX  operating 
system.  This  program  was  developed  to  the  point  of 
usefulness,  and  then  used  to  1  ic ate  events.  No  further 
program  development  was  done  (except  the  very  it.i  no'*  changes 


to  produce  f  ine  locate  and  far  locate )  .  so  there  are  sir* 


improvements  to  be  made. 

The  location  program  is  quite  interactive.  One 
hydrophone  with  a  bad  time  delay  can  change  the  slope  and 
location  a  great  deal.  An  event  is  located  by  eliminating 
bad  time  delays  and  checking  the  sound  speed  and  standard 
deviation  of  the  location.  In  some  cases  no  hydrophones 
needed  to  be  removed,  but  in  most  cases  at  least  one 
hydrophone  was  removed  before  an  event  was  considered 
located.  The  sound  speed  was  the  major  indicator  o-f 
whether  an  event  had  been  located.  If  the  sound  speed  was 
between  1380  and  1500  m/s  the  event  was  considered  located 
Of  course  an  attempt  was  made  to  get  close  to  1440  m/sec. 
This  had  to  be  balanced  with  reducing  the  standard 
deviation.  A  standard  deviation  below  0.01  seconds  was 
considered  good. 

A  table  summarising  all  of  the  events  and  their 
location  parameters  is  found  in  Appendix  C.  The  standard 
deviations  (sigma)  are  given  in  two  sets  of  units.  The 
first  is  the  sigma  calculated  by  the  location  program,  and 
it  is  in  seconds.  The  second  sigma  is  a  translation  of 
that  standard  deviation  to  meters  using  the  sound  speed 
calculated  for  each  particular  event.  The  standard 
deviations  ranged  from  0.0010  to  0.0327  sec,  with  0.0077 
sec  being  the  average.  The  significance  of  this  standard 
deviation  will  be  discussed  in  the  next  section. 

In  some  cases  just  removing  suspect  time  delays  did 


not  lead  to  a  1 ocal i sati on . 


A  reexamination  of  the  event 


time  series  was  done  to  see  if  any  of  the  manual  time 
delays  was  incorrect.  Often  a  reexamination  of  the  time 
series  produced  a  change  of  1  to  4  of  the  time  delays. 
These  corrected  values  plus  values  from  the  other  channels 
would  then  be  used  to  locate  the  event.  About  157.  of  the 
events  required  reexamination.  Most  of  those  were 
subsequently  located. 

Despite  the  above  efforts,  there  were  a  few  events 
that  could  not  be  located  within  the  1380  to  1500  m/s ec 
sound  speed  limits.  These  events  may  be  from  propagation 
paths  other  than  the  assumed  direct  acoustic  path.  Events 
arriving  primarily  through  the  ice  longitudinal  wave  or 
the  ice  flexural  wave  would  have  phase  speeds  above  and 
below  my  sound  speed  limits.  These  non-1 ocatabl e  evvents 
are  indicated  in  the  event  location  summary  of  Appendix  C, 
and  they  were  not  used  for  any  analysis  which  required 
accurate  location. 

Figure  3-5  shows  the  position  of  the  events  located 
within  a  2  km  square  centered  on  the  array  origin.  Figure 
3-6  shows  the  position  of  all  events  located. 
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Effects  of  Refraction  on  Location 


sound  speed  profile  was  used  to  get  the  refractive 


paths  for  various  ranges.  This  was  simplified  by  the  fact 


that  all  the  hydrophones  were  at  a  depth  of  93  meters. 


Assuming  that  only  the  "direct"  path  is  involved  means  that 


each  horizontal  range  has  to  have  a  unique  launch  angle  in 


order  to  reach  the  hydrophone  at  its  specific  depth.  Rav; 


were  1  a.i  rqhed  into  the  1  a  .  er 


the  sound  veiocit. 


and  the  horizontal  range  to  the  hydrophone  was  calculfi'ed. 


The  time  required  to  travel  the  refractive  path  can  be 


calculated  and  compared  to  that  of  the  slant  range.  This 


time  error  can  then  be  related  to  the  error  of  the  location 


program. 


Figure  3-7  shows  the  linearized  sound  velocity  profile 


that  was  used.  It  is  based  on  the  sound  velocity  orofile 


reported  for  the  eastern  Arctic  ocean  by  ChenCll.  Fiau 


-3  helps  to  illustrate  the  scheme  used  to  calculate  t 


ray  paths.  Equations  3.4, 


and  3.6  were  used  to 


calculate  angles,  ranges,  depths,  and  propagation  time. 


g  cos  B 


cos  Bq  -  cos  B  ^ 


•q  cos  B  _ 


iin  8^  -  sin  0. 
o  1 


t  =  —  I  1  n 


(1  •*-  sin  S^Ml  -  sin  Q^) 


(l  -  sin  0  .  )  (  l  +  sin  0  > 


Using  the  known  depths  (z)  and  estimated  speed 
gradients  (g)  of  Figure  3-7,  and  choosing  a  particular 
launch  angle  (0Q) ,  all  of  the  subsequent  angles  of 
intersection  of  the  layer  interfaces  (0^  ,  8^-,  0^)  can  be 
found  from  Equation  3-4.  The  angle  which  intercepts  the 
hydrophone  at  its  depth  (02>  can  also  be  found.  With  the 
angles  known,  the  horizontal  range  and  propagation  time  for 
each  layer  can  be  determined  with  Equations  3-5  and  3-b. 
These  are  combined  to  get  the  total  horizontal  range  and 
propagation  time.  It  should  be  noted  that  there  are  two 
ranges  and  times  for  each  launch  angle.  The  first  path  is 
that  which  intercepts  the  hydrophone  on  the  way  down,  while 
the  other  intercepts  the  hydrophone  as  it  is  refracted  back 
toward  the  surface.  The  maximum  depth  reached  by  the 
propagation  path  was  also  found,  and  those  paths  that  went 
below  754  meters,  resulting  in  a  range  greater  than  30600 
meters  were  not  reported.  A  tabular  summary  was  made  of 
launch  angle,  horizontal  range  from  the  hvdrophone,  maximum 
depth  and  propagation  time,  and  this  may  be  found  in 
Appendix  D,  along  with  mare  detailed  tables  listing  S^_^, 

r 1— 5  and  tl-4' 

Rays  connecting  source  and  hydrophone  with  one  or  more 
bounces  from  the  ice  were  not  considered  here.  The  effect 

of  those  rays  will  be  taler,  into  account  in  Chapter  . 

The  refractive  propagation  time  was  less  than  the 
slant  range  propagation  time  because  the  refracted  patn 


travels  through  faster  water.  The  slant  range  propagation 
time  was  calculated  by  dividing  the  slant  range  bv  1438.48 
m/sec,  the  average  sound  speed  between  0  and  93  meters 
depth.  The  time  di-f-ference  between  the  slant  range  path 
and  the  retracted  path  are  shown  in  Figure  3-9  as  a 
function  of  horizontal  range. 

This  time  difference  is  greater  than  the  average 
standard  deviation  of  the  location  program  or 1 j  after 
13,000  meters,  and  the  time  difference  at  20,000  meters  is 
only  about  twice  that  average.  The  standard  deviation  does 
not  reflect  the  time  difference  due  to  refraction  because 
all  of  the  hydrophone  time  delays  are  adjusted  in  the  same 
manner  and  direction.  Figure  3-10  shows  that  sigma  does 
not  grow  with  horizontal  range.  Refraction  e+fects  flc  not: 
influence  the  standard  deviation  greatly.  Closer  than 
13.000  meters  the  range  error  caused  bv  other  factors  masks 
any  error  from  ignoring  refraction. 

There  is  a  better  point  of  focus  for  examining  t re 
effect  of  refraction,  and  that  is  the  change  in  time  delay, 
not  the  change  in  the  propagation  time  itself.  A  point  was 
chosen  at  approximately  5  km  from  the  origin  of  the 
hydrophone  array,  and  another  chosen  at  approximately  6  km. 
The  slant  range  propagation  times  and  refractive 
prec-agst  i  on  *•  i  t-  es  were  calculated  for  each  point.  The  f  -.e 
delay  between  these  two  points  was  0.6857  sec  for  the 
•refractive  oath  and  0.6930  sec  for  the  slant  ranae  oath. 


Time  Difference  (sec) 


The  difference  was  0.0073  (about  1%)  ,  or  a  deviation  of 
0.0037  sec  -for  each  hydrophone.  Two  points  at 
approximately  IS  and  20  km  were  also  evaluated.  The 
di-rV  "~ence  between  their  time  delays  was  0.0182  sec  (about 
0.  V.)  ,  or  0.0091  sec  per  hydrophone.  These  numbers  are  the 
same  order  as  the  total  error  of  the  location  program. 

The  location  program  may  compensate  for  some  of  this 
error  by  raising  the  sound  speed.  If  just  the  points  above 
were  used,  the  sound  speed  would  go  from  1438.5  to  1453.8 
at  5  km  and  to  1458.6  at  20  km.  With  24  time  delays  being 
used  in  the  location  program  the  effect  may  not  be  as 
great. 

The  main  source  of  error  in  the  location  program  is 
the  quality  of  the  manual  time  delays.  When  the  signal-to- 
noise  ratio  was  low,  picking  the  correct  peak  was  often 
difficult.  The  standard  deviation  will  reflect  the 
judgement  of  the  person  picking  off  the  time  delays.  The 
time  delays  were  only  measured  to  the  closest  0.003  sec. 

It  is  interesting  to  note  that  0.006  seconds  equates  to  the 
width  of  a  pencil  tip  on  the  time  series  plot  scale. 

The  final  question  to  be  answered  is  "How  do  the 
standard  deviation  and  refraction  errors  equate  to  the 
range  and  bearing  accuracy  of  the  location  program. "  Two 
hypothetical  noise  events  were  investigated,  one  at  5000 
meters  (2845,  4136)  and  the  other  at  20,000  meters 
(-8253,  18896)  The  time  delays  for  slant  range  propagation 


and  refractive  propagation  were  calculated.  The  location 
program  was  run  -for  each  set  of  time  delays,  and  for  each 
set  partially  contaminated  with  0.016  sec  errors.  (Zero, 
+0.016  and  -0.016  were  each  added  to  one-third  of  the  time 
delays.)  The  results  are  summarized  in  Table  3—1. 


Table  3—1  Results  of  location  Program  Accuracy  Test 


Ar  o  A^  &  c 

(m)  (deg)  (deg)  (sec)  (m/s) 


slant  range 


5  km 


20  km 


ref ract i on 


20  km 


slant  range  w/ 


JO  km 


refraction  w/ 


5  km 


20  km 


4725 


18477 


4724 


18477 


4191 


20526 


4209 


20526 


-295 


55.6  +0.1  0.0003  1440 


-2143  113.7  +0.1  0.0002  1437 


-296  55.6  +0.1  0.0003  1454 


-2143  113.7  +0.1  0.0002  1458 


55.8  +0.3  0.0127  1453 


-94  113.9  +0.3  0.0130  1443 


-811  55.8  +0.3  0.0130  1466 


-94  113.9  +0.3  0.0130  1463 


The  refraction  contaminated  by  errors  case  is  closest 
to  what  was  input  into  the  location  program  for  the  field 


events , 


This  table  gives  an  estimate  of  the  accuracy  o f 


the  location  program  as  800  m  at  5  km,  and  2000  m  at  20  km. 

The  bearing  accuracy  is  excellent. 


i 


CHAPTER  4 


STRENGTH  OF  NOISE  EVENTS 

Acoustic  Source  Model 

The  dipole  is  considered  a  passible  source  model. 

Peak  values  Tor  the  source  parameter  of  -force,  F  ,  are 
used . 

The  acoustic  pressure  due  to  a  non— convect i ng  compact 

d  i  ode  source,  in  a  nonrefracting  infinite  medium,  isLSj; 


where  R  =  slant  range. 


Figure  4-1  shows  the  orientation  of  the  presumed  dipole. 

The  s.  gle  S  is  the  launch  angle  from  the  horizontal  plane 
down  into  the  water. 


the  denominator-  For  the  lowest  -frequency  considered  in 
this  study  <20  Hz)  ,  c  <1440  m/sec)  is  1 07.  of  2  77"  f  R  at  115 
meters  and  only  1.171  of  27TfR  at  1000  meters.  When  the  c 
in  the  denominator  is  neglected  the  force  can  be  written 


2p0R\ 


F  =  - 

°  sin  S 


(4-4i 


where  ^  is  the  wavelength. 


The  peak  pressure,  pQ  ,  should  lead  to  the  peak  force, 
Fq  .  This  definition  of  force  was  used  as  the  parameter 
for  dipole  strength.  Event  signatures  that  were  recorded 
from  a  source  withirr  300  m  of  a  hydrophone  were  not  usee  to 
calculate  dipole  strength,  FQ  ,  from  peak  pressure,  p0  . 

For  this  model  the  peak  acoustic  pressure  must  be 
found.  The  hvdroohcne  sensitivity  of  -15?  dB  re  1  volt  osr 


1  fj_  Fa  was  used  to  convert  voltage  to  pressureLlTl 


1  volt  =>  89  N/m^  =  89  Pa 


■-  4-5) 


The  dipole  strength  formula  requires  wavelength-  \ 
Frequency  was  taken  from  the  time  series  plots  for  each 
event  via  a;:  i  s  crossing  rate,  and  \  was  determined  bv 


dividing  c  <1440  m/sec)  bv  the  frenuenev. 

L. surer,  a.-.q  a  needed  for  ti  e  dioole  model  car,  p.  «•  c.,f 
as  in  Chapter  3  by  assuming  a  sound  velocity  profile  and 


imputing  the  refractive  path.  The1' 


unique?  1  amcl 
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angle  for  each  horizontal  range  when  the  path  is  purely 
retractive,  but  a  range  of  launch  angles  when  surface 
reflection  paths  are  included. 

The  final  parameter  in  the  dipole  strength  formula  is 
slant  range.  The  answers  obtained  using  slant  range  are 
the  strengths  based  on  spherical  spreading  in  a  nan- 
absorptive  medium,  equation  4-4.  Because  the  spherical 
spreading  assumption  is  a  poor  one,  ‘refractive  end  sur+a  :e 
reflective  propagation  paths  are  caused  by  the  Arctic  sound 
velocity  profile) ,  equation  4—4  must  be  modified.  The 
effect  of  refraction  on  spreading  loss  will  be  discussed  in 
the  next  section. 

Volumetric  absorption  was  found  by  using  the 
absorption  formulas  of  DyerC53-  Assuming  a  pH  of  8.2,  a 
salinity  of  33.5  °/oo,  a  temperature  of  0°  C  and  a  pressure 
of  40  atmospheres,  I  calculated  the  total  volumetric 
absorption  to  be  1.3  x  10  dB/km  for  an  80  Hz  signal.  For 
my  maximum  horizontal  range  of  20  km,  the  absorption  wculd 
be  0.026  dB.  This  is  not  significant,  and  I  therefore  did 
not  include  a  volumetric  absorption  correction  in  the 
strength  calculations. 

Effects  of  Refraction  on  Transmission  Loss 

Spherical  spreading  loss  in  a  ncnr e - r ac t : o a  medium  ; = 
illustrated  in  Figure  4-2C143.  The  sound  pressure  is* 
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The  sound  pressure  squared 
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is  proportional  to  intensity,  and  intensity  is  the  power 
per  unit  area.  Since  the  power  -from  a  source  is  const  an  r 

P  —  1  j  4  /  (  R  ^  =  I  *p4 7T  ^9^"  ■  ( 4-6  ) 

The  intensity  at  the  reference  range  o-f  1  meter  can  be 
related  to  other  intensities  by 


1 R 


^ref4^"  *ret 

47TR2  R2 


,  4  —  7  ) 


Since  I  =  p'/pc  ,  this  can  be  expressed  in  terms  o-f 
transmission  loss,  H  . 


n 

PR^  ^ 

H  =  -10  log  - —  =  10  log  R^  =  20  1 og  R  ,  (4-B) 

°r  e-f 

in  dB  re  the  distance  reference,  taken  as  1  m. 

The  spreading  scheme  tor  a  retractive  medium  is  shown 
in  Figure  4-3C23.  This  is  based  on  ray  theory  which 
assumes  that  acoustic  energy  does  not  cross  the  ravs,  with 
energy  contained  between  two  rays  being  conserved.  The 
intensity  at  the  reterence  range  between  the  two  rays  shown 
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P 


2 77" R c 1  j s  9d  R^s  ZTTcos  bq 


(4-9  ) 


At  a  horizontal  distance  r  meters  from  the  source,  the 


l n tens i t y  is: 


<4-1 0 ) 


I,  = 


2  77>L  2JT  r  /Sr  51  n  aTT  27T r  ^rcos  81 


The  relation  between  intensities  becomes: 


Ire+27Tcos  8a^3  _  Ire^^Bcos  SQ 

r  27T r  ^rl5in  Sll  r^rjsin  Sj 


(4-1  1 ) 


The  loss  due  to  spreading  is: 


_  Ae 


Qcos  a. 


rret  rAr,5in  S ^ 


(4-12) 


In  terms  of  pressure,  Ir  =  p^/y^^c^,  and  Ire^ 
and  there-fore 


pre+  /  Poc 


Pici  A0  c°5  dr- 


ret 


Poco  r  Ar|Sin  81 


4-  1 


hi  nee  b v  Snell's  law.  cos  0_  ■  c_,  =  cos  9.  c. 

o  o  i 


.  Pi  A 


S  cos  9, 


A« 


pref 


r  £r  |sin  Sj  r^rjtan  sj 


since  /  yQ  Q  =  1  in  seawater  tr  an 

excel  1  ent  apor  o  .■  me*-  ;  -r, 


Appl  ying  the  dipole  model  to  this  spreading  l 
equation  qives 
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Equation  4—15  assumes  a  unique  refractive  path  between 


source  and  hydrophone.  In  the  Arctic  there  may  be  other 
paths  due  to  the  non— specular  scattering  of  rays  off  the 
ice  canopy.  Figure  4-4  illustrates  how  rays  normally 
trapped  in  a  surface  duct,  may  be  deflected  down  to  a 
hydrophone.  The  minimum  vertexing  angle  calculated  from 
the  linearized  sound  velocity  profile  of  Chapter  3  was 
0.064  radians.  For  a  ray  to  stay  in  a  surface  duct  above 
the  hydrophone  it  must  be  reflected  from  a  slope  of  less 


than  0.032  radians  or  about 


It  is  reasonable  to  assume 


that  the  ice  canopy  lacks  local  levelness  to  this  order,  so 
that  non— specular  rays  must  be  accounted  for. 

The  rays  which  rebound  from  the  ice  canopy  experience 
some  loss.  The  attenuation  for  the  FRAM  IV  experiment  has 
been  reported  at  0.1  dB/km  at  80  Hz  C 1 1 ]  .  This  attentuation 
may  be  converted  to  a  loss  per  bounce. 


P  * 


0. 1  dB/km  = 


(4-1  to) 


where  b  =  loss  per  bounce,  and 
X  =  cycle  distance. 


The  cycle  distance  depends  on  the  launch  angle  and  the 
sound  speed  gradient.  For  a  launch  angle  of  0.032  radians 
and  the  assumed  sound  velocitv  profile  of  Chaoter  3.  the 
cycle  distance  is  3.7  km.  Therefore,  the  loss  per  bounce 
is  about  0.4  dB.  This  loss  is  low  enough  that  even  a  ray 
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which  has  bounced  several  times  may  contribute  a 
significant  amount  of  energy  at  the  hydrophone.  The  non- 
specular  rays  cannot  be  ignored.  Bounce  loss  at 
frequencies  less  than  80  Hz  are  even  smaller,  since  the 
data  show  a  roughly  linear  dependence  on  frequency. 

To  account  for  the  non-specul ar  rays  the  spreading 
loss  is  calculated  using  the  ray  averaging  techniqueCSl . 
The-  pressure  from  a  particular  ray  at  a  given  depth  and 
horizontal  range,  assuming  a  dipole  source  model,  is 


p~(r,z)  = 


51n^8o  d8  dr 
r | tan  |  dr  X/2 


(4-17) 


The  term  -  represents  the  probabilitv  that  a  rav  bundle 

X/2 

will  cross  a  certain  depth,  as  shown  in  Figure  4-5.  For  a 
single  linear  sound  speed  gradient  the  cvcle  distance  can 
be  written  as 


sin  eo  = 


2r_  sin  9„ 
c  o 


2rcSo 


'.4-18) 


where  rc  is  the  radius  of 
curvature,  to  a  good  approximation  constant 
for  all  small  angle  rays  in  a  linear  sound 
speed  gradient. 


Applying  equation  4-18  to  equation  4-17,  and  using  the 
small  angle  aopr  o>:  i  mat  i  on  .  gives 


P*~  <r  ,z  )  = 


*2  leol  d8o 


(4-19) 


rc  r 


f •.*< 


Figure  4-5  Probability  of  a  ray  bundle  crossing  a  certain 
depth  at  a  given  horizontal  range. 


vs- v>.-v- 

^  ^  1  _  •  .  v_  4  .  ■  _  ^  ,*  _ 


-68- 


In  order  to  average  the  contributions  of  the  possible 
rays,  this  pressure  is  integrated  over  all  possible  angles 
•for  a  given  receiver  depth,  and  then  averaged  over  depth 
down  to  the  hydrophone  at  z Q. 


where  @v  is  the  maximum  launch 
angle  o-f  a  ray  that  will  hit  the  hydrophone 
at  a  given  range,  and  Sm  is  the  minimum 
launch  angle. 


The  angle  is  a  function  o-f  8Q  and  z 

z  =  rc  Ceos  0^^  -  cos  B0D  .  (4-21 

Using  the  small  angle  approx  i  mat  i  on  -for  cosine  leads  to 


t  4-22 


Substituting  this  into  equation  4-20  and  evaluating  the 

integral  over  angle  gives 


0_ "  =  —  for  all  z  so  the  second  term  within  the  intecral 
rc 

is  always  zero.  Evaluating  the  first  term  over  deoth  ci  /e- 

an  expression  for  pressure  in  terms  of  r  and  B  . 
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This  expression  can  be  used  to  -find  the  source  strength. 


The  spreading  function,  G,  is  presented  as  a  function  of  r 
alone  since  Sv  depends  on  r.  For  each  r  there  is  a  unique 
9V ,  and  therefore,  a  unique  G.  The  spreading  function  was 
calculated  for  horizontal  ranges  from  300  m  to  20,000  m, 
and  tablulated  in  Appendix  D.  The  spreading  function  is 
shown  in  a  log-log  plot  in  Figure  4-6.  For  comDarison  the 
equivalent  spherical  spreading  for  a  dipole  source  is  also 
shown.  From  this  one  can  see  that  source  strengths 
calculated  using  the  spherical  spreading  law  lead  to  an 
unrealistic  dependence  on  range. 

In  order  to  get  8X/  and  rc  a  linear  sound  speed 
gradient  of  0.054  sec  1  was  chosen.  This  gradient  gives 
the  same  0S/  at  r  =  3  km  as  the  multiple  step  profile  used 
in  Chapter  3.  Three  kilometers  was  chosen  si.nce  it  was  the 
median  horizontal  range  for  the  noise  events. 

The  spreading  loss  function  and  measured  peak  pressure 


magnitudes  were  used  to  calculate  dinole  strength. 


CdB 


For  a  part  1  :ul  nr  event  'he  sr  ;rce  strength  was 
calculated  f  r  'im  ea<  h  u  ,  lr  I'inrip,  and  fhen  the  average  taten 
as  the  source  e 1  -  •  ’  ■  •  ••  •  •  ■<-  I  he  standard 

devi  at  jn  within  .=  -«  •  .  »-•  •  ■  «  *  -  >m  to  1 75'/.  at  the  mean 

V  a  1  J  .  no  '  ’  *  '  n  (  i  ( i  n  f  t  d  1  C  lists  the 

mean  measnr  e<i  e  . .  .  •  ■•*  w.m  dipole  strength, 

along  wi*h  •  •.«=.  ■  .  •.  /allies  tor  each 

,  ,  l 

Gt^-enijth  h,.  i  j'  •  • 

I  was  inte<-e<-’  i  ■  ■  •  *  Mat  =»n  /  i  r  on  men  t  a  1 

1  oad  inq  might  iij,p  had  in  *  ►  e  t  e»mpi  ir  a  1  ,  spatial  and 
s  t  r  eng  t  h  statistic.  ;  t  t.as  heen  V  •  vwn  h  /  Malris  and 
Dyer  CIO]  tnat  low  tregneni  /  1  1  Jo  Hz  band)  ambient  noise 
rms  pressure,  aver  aged  over  a  long  time,  correlates  well 
with  environmental  stresses  and  Tinmen'  s.  Since  I  had 
ambient  noi  ae  pres  sore  c  >r  must  ot  the  period  of  r.he  FRAId 
IV  experiment,  and  since  I  had  en v i ronmental  stresses  and 
moments  available  for  only  a  part  of  the  time,  I  chose  to 
use  the  20-80  Hz  long-time-average  rms  pressure  as  my 
envi ronmental  indicator. 

The  10  to  20  Hz  band  ambient  noise  pressure  was 
converted  to  20  to  80  Hz  band  pressure  in  the  following 
manner .  Figure  4-7  shows  the  typical  spectrum  for  central 
Arctic  pack  ice  noise.  The  portion  of  the  spectrum  between 
10  and  100  Hz  can  be  approx i mated  by  a  straight  line. 


At 


log  S  =  A  [  log  F  J  +  B  , 


(4-27 


where  A  =  slope 


=  -1.7273  Pa^/Hz- 


B  =  intercept  =  —1.0909  Pa^/Hz  , 


S  =  10BFA 


(4-28 


The  band  rms  oressure  relates  to  the  soectral  level  ov: 


rms ,  b 


S  dF 


(4-29 


I  have  assumed  that  as  the  sound  pressure  level  changes 
From  time  to  time  the  intercept  B  changes,  but  the  since 
remains  the  same.  By  substituting  equation  4-28  into  4-29 
and  using  the  known  10  to  20  Hz  ambient  noise  band,  B  can 
be  written  in  terms  of  the  known  oressure. 


B  =  1  og 


rms . 10—20 


(4-30 


where 


K'l  " 


9Q ^  ^  +  1 )  _  10(A 

A  +  1 


=  0.1 020  H: 


The  ambient  noise  rms  pressure  For  the  20  tc  30  Hz  band 
now  be  Found. 


rms . 20-80 


O  =  / 


Table  4-1 

20  to  SO  Hz  Band  Ambient 

Noise  rms  Pressure 

T ape  # 

Date  Recorded 

prms , 20-80 
(Pa) 

400 1 

3—27—82 

Not  Aval  1 ab 1 e 

2001 

3-29-82 

Not  Avai 1 abl e 

2009 

3-30-82 

0.022 

300 1 

3-31-82 

0.019 

4003 

4-01-82 

0. 044 

4005 

4-01-82 

0 . 035 

4007 

4-01-82 

0 . 022 

4009 

4-02-82 

0.01 0 

4011 

4-02-82 

0.010 

40 1 3 

4-03-82 

0.013 

2023 

4-08-82 

0.037 

4015 

4-09-82 

0.  040 

3047 

4-13-82 

0.010 

4016 

4-15-82 

0.017 

40 1 9 

4-15-82 

0 .016 

4021 

4-19-82 

0 . 0 1  3 

4023 

4-19-82 

0 .011 

4024 

4-19-82 

0 .011 

4027 

4-20-82 

0 .012 

4029 

4-20-82 

0 . 0  1  2 

4031 

4-20-82 

0 .012 

4033 

4-20-82 

0.012 

4040 

4-21-82 

0.  034 

4047 

4-21-82 

0.  114 

4049 

4-21-82 

0 .  1  40 

405  1 

4-21-82 

0.  140 

4053 

4-22-82 

0 . 080 

4055 

4_?n_Q2 

0 . 082 

4057 

4-22-82 

0. 053 

4059 

4 _ 79 _ gn 

0.065 

4061 

4 _ n n _ go 

0.034 

4063 

4—22—82 

0.0  23 

CHAPTER  5 


ANALYSIS  OF  NOISE  EVENTS 
Detection  Analysis 

A  total  of  34  tapes  was  examined,  -for  a  total  time  of 
662  minutes.  (For  a  few  of  these  tapes  the  entire  20 
minutes  was  not  used.) 

There  was  a  total  of  499  detections  of  events  flagged 
on  at  least  50'/.  of  the  hydrophone  channels.  Of  these,  139 
were  man-made  artifacts,  and  125  were  false  alarms 
(detections  which  were  so  weak  that  no  pattern  for  taking 
time  delays  could  be  discerned).  There  were  199  unique 
events,  and  36  multiple  dtections  of  those  events.  Stated 
in  another  way,  of  the  detections  which  were  not  artifacts, 
65.3%  were  strong  enough  to  support  analysis  and  34.7%  were 
too  weak  to  reasonably  analyze,  and  hence  labeled  false 
al arms. 

Since  the  detection  process  depends  on  si gnal -to-noi se 
ratio,  the  level  of  background  ambient  noise  should  affect 
the  event  detection  rate.  Figure  5-1  shows  normalized 
ambient  noise  pressure,  number  of  false  alarms  per  tape, 
and  number  of  unique  events  per  tape  for  each  tape 
examined.  There  is  some  trend  for  more  events  being  found 
when  the  ambient  pressure  is  low,  and  more  false  alarms 
■ declared  when  the  ambient  pressure  is  high. 

This  is  more  clearly  seen  in  Figure  5-2,  which  shows 
the  average  number  of  false  alarms  and  unique  events  found 
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per  tape  in  each  of  four  background  noise  pressure  ranges. 
The  0.01-0.02  Pa  range  used  15  tapes  to  compute  its 
average,  the  0.02-0.03  Pa  range  4  tapes,  the  0.03-0.04  Pa 
range  5  tapes,  and  the  over  0.04  Pa  range  8  tapes.  Two 
tapes  were  recorded  during  the  first  few  days  of  the  FRAU 
IV  experiemnt,  before  the  10-20  Nz  band  ambient  noise 
recordings  were  started. 

A  breakdown  of  detections  for  each  tape  is  found  in 
Appendix  C. 


Temporal  Analysis 

The  interarrival  time  between  events  ranged  from  1  to 
1064  seconds.  Each  event  time  was  taken  to  the  nearest 
second,  and  no  events  were  taken  as  having  the  same  event 
time.  If  two  events  happened  in  the  same  second,  one  was 
judged  to  be  earlier,  and  the  two  events  were  given  event 
times  one  second  apart.  The  interarrival  time  for  a 
particular  event  was  measured  from  the  previous  event, 
except  for  the  first  event  of  a  tape,  which  was  measured 
from  the  start  of  the  tape. 

The  interarrival  times  were  divided  into  bins  of  20 
seconds.  The  first  bin  ("0" )  contained  events  which  had 
interarrival  times  from  0  to  19  seconds,  the  second  bin 
from  20  to  39  seconds,  and  so  on.  The  number  of  events  per 
bin  is  presented  in  Table  5-1  and  shown  graphically  in 
Figure  5-3.  A  complete  listing  of  interarrival  times  for 


Number  of  Events 
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each  event 

is  -found  in  Appendix 

C.  The 

mean  o-f  the 

interarrival  times  ( [Ji )  is  10O 

seconds , 

and  the  standard 

devi ation 

<(J)  166  seconds.  In 

terms  o-f 

bins,  the  mean 

5  and  the 

standard  deviation  8. 

The  standard  deviation 

1.66  times 

the  mean. 

Table  5-1 

Number  o-f  Events  per 

Interarrival  Time  Bin 

Bi  n 

Events 

Bin 

Events 

0 

67 

1 

36 

31 

0 

O 

<Lm 

19 

32 

0 

3 

15 

33 

1 

4 

12 

34 

0 

5 

7 

35 

0 

6 

4 

36 

0 

7 

3 

37 

0 

8 

6 

38 

o 

9 

1 

39 

0 

io 

3 

40 

0 

1  1 

4 

41 

1 

12 

2 

42 

0 

13 

1 

43 

0 

14 

3 

44 

0 

15 

2 

45 

0 

16 

0 

46 

0 

17 

0 

47 

0 

IB 

0 

48 

0 

19 

0 

49 

1 

20 

1 

50 

0 

21 

0 

51 

0 

22 

0 

52 

0 

23 

2 

53 

1 

24 

1 

54 

0 

25 

1 

55 

0 

26 

1 

56 

0 

27 

2 

57 

0 

28 

1 

58 

0 

29 

1 

59 

0 

30 

0 
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Three  different  probability  density  functions  were 
investigated  to  find  an  appropriate  fit  for  Figure  5-3. 
They  were  1)  a  hal f -gaussi an  distribution,  2)  an 
exponential  distribution  and  3)  a  J  shaped  distribution. 

The  hal f -gaussi an  probability  density  function  isC31: 


p  (t ) 


<5-1 ) 


The  general  equations  for  mean,  mean  square  value  and 
variance  ((72)  can  be  used  to  solve  for  the  unkown 
constant,  tQs 


■/ 


CO 


t  p (t )  dt 


(5-2) 


mean  square  value  -  J  t ^  p(t)  dt 

0 


-/ 


CO 


(5-3) 


•/ 


CO 


(J  ‘**  =  J  (t  ~  fJL  >  ^  p<t)  dt 
0 


(5-4) 


=  mean 


square  value  - 


Substituting  equation  5-1  into  equations  5-2,  5-3  and  5-4 
leads  to  the  following  relations: 


y 

mean  square  value  =  tQ 
tQ2  ;  (J=  0.756  jJi 


(5-5) 


This  value  -for  tQ  was  used  in  equation  5—1,  and  the 
probability  density  -function  integrated  over  appropriate 
limits  to  get  the  number  o-f  events  in  each  20  second  bin. 
The  result  is  plotted  against  the  experimental  distribution 
in  Figure  5-4. 

The  second  distribution  (the  exponential)  belongs  to 
the  family  of  gamma  distribution  functions  C153: 


P 


1 

to  O' +  1  r<  a  + 


a  -t/tQ 

t  e 


(5-6) 


When  =0,  this  becomes  the  exponential  probability 
density  function 


P  = 


-t/t0 


(5-7) 


Again  using  equations  5-2,  5-3  and  5-4  leads  to: 


2 

fj.  =  tQ  ;  mean  square  value  =  2 fj, 

cr2  -  ijl2  i  <7 -ii. 


(5-8) 


The  exponential  probability  density  function  was  integrated 
over  the  bins,  and  the  results  are  shown  in  Figure  5-5. 

Another  demonstration  of  the  fit  of  the  exponential 
probabilty  distribution  is  shown  in  Figure  5-6.  Taking  the 
natural  log  of  the  function  should  lead  to  a  straight  line 


when  plotted  against  time  or  bin  number. 


The  straight  line 
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in  Figure  5-6  is  a  plot  of  the  natural  log  of  the  points 
calculated  using  the  exponential  probability  density 
•function-  The  experimental  points  seem  to  curve  rather 
than  lie  on  a  straight  line. 

The  last  distribution  <J  shaped)  is  also  a  gamma 
distribution.  The  J  shaped  distributions  are  characteri zed 
by  a  <  0  .  I  chose  a  -fairly  common  distribution  with 
Qf  =  -0.5  .  The  probability  density  function  is: 

i  -1/2  -t/tQ 

p(t)  =  -  t  e  ,  (5-9) 

V^T 

and  the  key  parameters  are: 


mean  square  value 


<7-V^U  ■ 


( 5- 1 0 ) 


This  distribution  is  plotted  against  the  experimental 
values  in  Figure  5-7.  The  natural  log  of  both  calculated 
and  experimental  points  are  plotted  against  bin  number  in 
Figure  5-8.  This  distribution  seems  to  fit  the 
experimental  points  best  of  all.  The  J  shaped  probabi ltv 
density  function  goes  to  infinity  at  zero,  but  it  is 
in'-.egrable. 

A  Chi  square  goodness  of  fit  test  was  done  on  all 


three  distributions. 


The  results  are  summarized  in  Table 


Number  of  Events 


0  2  4  6  8  10  12  14  16 

20  Second  Bin  (0=0- 19s,  1=20-39s,  etc) 


Figure  5-7  J  shaped  distribution  compared  with 
experimental  values. 


(number  of  events) 
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5—2.  Also  presented  in  Table  5—2  are  the  ratios  of 
standard  deviation  to  mean. 

Table  5—2  Comparison  o-f  Distribution  Functions 


Chi  square 

( J//JL 

Ex per i mental 

— 

1.66 

Hal f — gaussi an 

127.  17 

0.76 

Exponenti al 

45.69 

1 . 00 

J  shaped 

10.65 

1 .41 

For  a  distribution  to  pass  a  goodness  o-f  -fit  test  it 
must  have  a  Chi  square  less  than  a  prescribed  limit.  The 
limit  -for  my  test  (9  degress  o-f  -freedom,  Ql  =  0.005)  was 
23.6C161.  Only  the  J  shaped  distribution  passed  the  Chi 
square  test.  It  also  has  (J  /  U  closest  to  the 
experimental  values.  In  summary,  the  interarrival  data 
reasonably  fit  a  J  shaped  distribution  given  by: 

1  -1/2  -t/2 U 

p(t)  =  -  t  e  .  (5-11) 

Since  event  detection  rate  depended  on  ambient  noise 
level,  interarrival  time  between  events  should  also  show 
environmental  dependence.  Table  5.3  gives  average  and 
standard  deviation  of  the  interarrival  time  for  different 
ambient  noise  pressure  levels. 
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Table  5-3  Background  Noise  Level  Dependence  o-f 

Interarrival  Time 


Ambient  Noise  Interarrival  Time 

rms  pressure 


(20—80  Hz) 

(Pa) 

mean 

( sec ) 

standard  deviation 

(sec) 

0.01-0.02 

67 

129 

0.02-0. 03 

190 

153 

0.03-0.04 

147 

143 

over  0.04 

183 

318 

The  tapes  having  a  background  noise  level  of  0.01  to  0.02 
Pa  have  a  significantly  shorter  inter arrival  time  than 
tapes  in  the  other  three  pressure  groups.  As  with 
detection  rate,  the  inter arr ival  tine  does  depend  on 
ambient  noise  level. 

Spatial  Analysis 

After  removing  nonlocatable  events  and  events  located 
outside  a  horizontal  range  of  20,000  meters,  164  events 
remained.  These  were  grouped  by  horizontal  range  into  42 
annuli  of  equal  area  as  shown  in  Figure  5-9.  Each  annulus 
is  a  30  square  km  ring  centered  at  the  array  origin.  The 
first  annulus  ("0")  went  from  0  to  3090  meters,  the  second 
from  3090  to  4370  meters,  and  so  on. 

Table  5—4  shows  the  number  of  events  per  annulus  and 
Figure  5-10  shows  this  distribution  graphically. 


Table  5—4  Number  at  Events  per  Annulus 


inulus 

Events 

Annul  us 

Even' 

0 

91 

21 

0 

1 

19 

n  n 

0 

ni 

4 

n  ~r 
•-> 

1 

,j 

24 

0 

4 

O 

4. 

25 

0 

5 

4 

26 

1 

6 

9 

27 

0 

7 

3 

28 

0 

a 

1 

29 

0 

9 

cr 

30 

o 

10 

n 

31 

Cj 

1 1 

3 

*T 

0 

12 

4 

1 

13 

1 

34 

0 

14 

1 

35 

1 

15 

4 

36 

0 

16 

0 

37 

1 

17 

0 

38 

1 

18 

1 

39 

0 

19 

0 

40 

0 

20 

o 

41 

1 

The  average  number  at  events  per  annulus  is  3.93  and  the 
standard  deviation  is  14.15  events.  Figure  5-10  shows  that 
the  number  of  events  found  is  highly  dependent  on  their 
range  from  the  array.  In  the  center  annulus  there  were 
over  20  times  the  mean  number  o-f  events. 

The  dependence  on  range  is  not  a  surprise,  since 
spreading  (and  possibly  scattering  and  other  losses)  will 
reduce  the  strength  at  weal  transients  down  to  the  ambient 
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square  kilometer  per  hour  . 

The  average  number  o- f  events  per  annulus  and  the 
number  of  events  per  square  kilometer  per  hour  should 
depend  on  background  noise  level.  The  average  number  of 
events  per  annulus  was  -found  -for  each  ambient  noise  rms 
pressure  range,  and  adjusted  to  reflect  the  number  of 
events  in  a  662  minute  period.  The  results  are  seen  in 
Table  5—5,  The  number  of  events  per  square  kilometer  per 
hour  for  the  center  annulus  are  also  shown  in  Table  5-5. 

Table  5-5  Average  Number  of  Events  per  Annulus 


for  4  Ambient  Noise 

Level  s 

Ambient  Noise 
rms  pressure 

(20-80  He) 

Events 

per 

Annul  us 

Mi nutes 
of  tape 
Ex  ami  ned 

Ad  justed 
Events  per 

Annul  us 

#  Event 
per  km 
per  hr 

0.01-0.02  Pa 

2.60 

287.5 

5.98 

0. 452 

0.02-0.03  Pa 

0.29 

77 

2.46 

0.  260 

0.03— O. 04  Pa 

0.  45 

1  00 

2 . 99 

0.  130 

over  0.04  Pa 

0  -  36 

157.5 

1 . 50 

0. 076 

Ent  l  re 

Population  3.93  662  3.93  0.275 

The  average  number  of  events  per  annulus  and  the 
number  of  events  per  square  kilometer  per  hour  bi>t  h  reflect 
the  effect  of  s i qn a  1 -to-no i se  rati  >  on  the  detection 

-  ■  :'i •  - me  . 

The  entire  papulation  of  events  was  investigated  for 
angular  dependence.  Figures  5-11  and  5-12  show  the  number 
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Figure  5-11  Number  o-f  events  per  30°  sector.  Ancles  are 
measured  -from  the  northern  leg  o-f  the  hvdrochone  arrav. 


of  events  -found  per  30  sector.  In  the  polar  diagram 
(Figure  5-12)  the  radius  shows  the  number  o-f  events.  The 
angles  are  measured  -from  the  northern  leg  o-f  the  array. 
Figure  5—13  is  a  polar  plot  showing  the  number  of  events 
per  10°  sector.  There  was  no  predominant  angular  direction 
found.  However,  some  preference  can  be  seen  -for  bearings 
o-f  330°  and  190°  -from  the  northern  leg  o-f  the  array. 


Strength  Analysis 

The  mean  hydrophone  peak  pressure  magnitude  for  each 
event  fell  within  a  fairly  narrow  band  of  values.  The  mean 
peak  pressures  ranged  from  1.32  to  0.16  Pa,  with  an  average 
of  0.36  Pa  and  a  standard  deviation  of  0.20  Pa.  Figure  5-1 
shows  the  mean  hydrophone  peak  pressure  values  for  all 
events  located  between  100  m  and  20,000  m  plotted  against 
range  from  the  array  origin. 

The  different  symbols  shown  in  Figure  5-14  represent 
events  during  each  of  the  four  ambient  pressure  categories. 
The  events  with  a  higher  mean  hydrophone  peak  pressure  have 
a  tendency  to  occur  during  higher  ambient  noise  levels. 

This  can  be  seen  in  Table  5-6,  where  the  maximum,  minimum, 
average  and  standard  deviation  of  the  mean  hydrophone  peak 
pressure  values  are  given  for  each  of  the  four  ambient 
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Figure  5-14  Mean  hydrophone  peak  pressure  measured  for 
events  between  100m  and  20,000  m,  plotted  against 
horizontal  range  from  the  FRAM  IV  array  origin. 
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Table  5—6  Mean  Hydrophone  Peak  Pressure 
■for  4  Ambient  Noise  Levels 


Ambient  Noise 

Mean 

Hydrophone  Peak  Pressure 

rms  pressure 
( 20-80  Hz ) 

max 

(Pa) 

mi  n 

average 

std  dev 

0.01-0.02  Pa 

1.32 

0.  16 

0.31 

0.17 

0.02-0.03  Pa 

0 . 56 

0.  25 

0.  36 

0.09 

0.03-0.04  Pa 

0.85 

0.  23 

0.  43 

0.  15 

over  0.04  Pa 

1 . 28 

0.  16 

0. 49 

0 . 34 

Ent i re 

Popul at i on 

1 . 32 

0.  16 

0 . 36 

0 .  20 

Source  strength  (FQ) 

was  found  for 

the  events 

whi  ch 

had  hydrophone 

locations 

between  300  m 

and  20,000  m 

f  r  om 

the  event. 

The  dipole  strengths  ranged  from  33  kN  to  4.9  MN , 
with  an  average  of  431  kN  and  a  standard  deviation  of  555 
kN.  The  distribution  of  strengths  for  the  151  events 
evaluated  is  shown  in  Figure  5-15  and  in  Table  5-7. 

Figure  5-16  shows  the  dipole  strength  for  all  events 
plotted  against  horizontal  range  from  the  array  origin. 
Again,  it  can  be  seen  that  the  stronger  events  occur  w hen 
the  ambient  pressure  level  is  high.  Table  5-8  oives  the 
strength  values  for  the  different  ambient  noise  levels. 


Table  5-7  Strength 

Distribution  -for 

a  Population 

of  Events 

Fo 

#  o-f  Events 

0 

to 

100 

kN 

19 

100 

to 

200 

kN 

29 

200 

to 

300 

kN 

33 

300 

to 

400 

kN 

28 

400 

to 

500 

kN 

8 

500 

to 

600 

kN 

5 

600 

to 

700 

kN 

4 

700 

to 

800 

kN 

5 

800 

to 

900 

kN 

2 

900 

to 

1000 

kN 

5 

1000 

to 

1100 

kN 

3 

110O 

to 

1200 

kN 

2 

1200 

to 

1300 

kN 

3 

1300 

to 

1400 

kN 

1 

1400 

to 

1500 

kN 

1 

1.5 

to 

3 

MN 

1 

over 

3 

MN 

2 

Table  5-8 

Di pal e 

Strength 

versus 

Ambient  Noise 

Level s 

Amb i en t  Noi se 

Di pole 

Strength  (FQ) 

rms  pressure 

(kN) 

(20-80  Hz) 

max 

mi  n 

average 

std  dev 

0.01-0.02  Pa 

1051 

33 

259 

155 

0.02-0.03  Pa 

2041 

64 

649 

634 

0.03-0.04  Pa 

1 153 

79 

643 

359 

over  0.04  Pa 

49 

’39 

59 

860 

1449 

Entl re 

Popu 1  at i on 

49 

’39 

_«.j 

431 

555 

Horizontal  Range  From  Array  Center  (m) 
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igure  5-16  Dipale  strength  far  events  between  300  m  and 
0,000  m,  platted  against  horizontal  range  from  the  array 
ir  i  gi  n . 
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Most  of  the  events  evaluated  -for  strength  occurred 
during  the  lowest  ambient  noise  levels.  Figure  5—17  shows 
the  strength  of  events  that  occurred  when  the  ambient  noise 
was  0.01  to  0.02  Pa.  The  log  o-f  the  dipole  strength  is 
plotted  against  the  log  of  the  horizontal  range  from  the 
center  of  the  array.  The  points  scatter  more  so  to  the 
upper  left  rather  than  lower  right,  because  distance  itself 
filters  out  weak  events.  A  weak  signal  from  far  away  would 
not  reach  the  hydrophone  array  with  enough  amplitude  to  be 
distinguished  from  the  background  noise.  And  events 
located  farther  away  would  tend  to  be  strong  events. 
However,  events  located  close  to  the  array  should  have  the 
entire  range  of  source  strength  levels.  This  would  produce 
a  wedged  shaped  plat  of  weaker  events  close  to  the  array. 
Indeed,  Figure  5-17  shows  a  general  scattering  with  perhaos 
a  wedge  of  weaker  events  near  the  array  origin. 

Nonetheless  the  trend  shown  in  Figure  5-17  suggests 
that  the  ray  average  model  used  to  estimate  refractive- 
surface  reflective  spreading  may  need  to  be  replaced  with  a 
more  refined  model.  For  example,  horizontal  ranges  less 
than  about  1000  m  may  include  too  smal.l  a  loss,  and 
therefore  lead  to  too  small  a  strength,  because  the 
reflective  contributions  mav  not  be  as  large  as  imputed. 
Such  a  criticism  is  supported  b y  the  notion  that  tc-r  a 
given  slope,  £  ,  reflective  rays  and  hence  rav  averaging 

occurs  only  beyond  a  critical  horizontal  ranqe. 
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The  foregoing  speculation  suggests  that  the  average 
dipole  strength  for  the  lowest  ambient  noise  case  is  best 
found  from  the  events  farther  from  the  array,  and  is 


F  *  103-°  N  *  320  kN  , 


(5-12) 


with  a  much  smaller  standard  deviation  than  in  Table  5-8. 
Presumably  correspondi ng  adjustments  could  be  made  for  the 
higher  ambient  noise  cases,  but  the  FRAM  IV  data  set 
contains  too  few  events  at  higher  ambient  noise  to  plot  as 
in  Figure  5-17. 

The  strength  analysis  is  a  somewhat  ambivalent  one 
because  of  spreading  model  uncertainty,  and  because  data  on 
ice  slopes  are  not  available.  But  the  dipole  picture  of  an 
event  likely  has  some  validity,  and  at  least  rough 
estimates  of  its  strength  have  been  extracted  from  the 
data . 
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CHAPTER  6 

SUMMARY  AND  THOUGHTS 

Through  the  use  of  a  detection  program,  visual 
confirmation  and  a  location  program,  a  population  of  199 
Arctic  noise  transients  was  gathered.  There  are  -four  major 
resul ts. 

First,  more  events  are  found  when  the  ambient  or  ess  ire 

is  low,  and  more  false  alarms  when  the  ambient  pressure  is 
high.  The  interarrival  time  and  the  average  number  of 
events  per  unit  area  also  depend  on  ambient  noise  level. 
Since  mare  event-s  are  found  when  the  ambient  noise  is  low, 
the  interarrival  time  decreases,  and  the  spatial  density 
i ncreases. 

Second,  the  inter arrival  times  were  fit  to  several 
possible  probability  distributions.  The  interarrivai  time 
distribution  best  fits  a  J  shaped  gamma  distribution.  The 
mean  interarrivai  time  is  100  seconds. 

Third,  the  number  of  events  per  unit  area  is  highly 
dependent  on  range,  since  distance  filters  out  weak 
transients.  The  event  density  in  the  annulus  closest  to 


the 

center  of 

the 

array  was  0.3  events 

per  square 

k  i  1  ofneter 

per 

hour  over 

al  1 

observations  and  0.5 

events  per 

saucer  e 

i  ■» 

*  1  A  l 

meter  e  r 

hour 

for  aui et  times.  ~h 

ere  l  =  r.o 

predominant  angular  dependence  to  the  spatial  distribution 
of  events. 
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Last ,  the  mean  dipole  strength  tor  the  observed  events 
is  430  kN  overall  and  260  1 N  during  low  ambient  noise 
levels.  Stranger  events  occurred  during  high  ambient  noise 
levels.  A  retinement  of  the  spreading  loss  model  used  to 
calculate  these  values  may  lead  to  values  which  are 
slightly  higher. 


Hiiaiyais  ot  hi  ccic  acoustic  events  is  tar  trom 

complete.  Several  areas  tor  improvement  have  been 
mentioned  earlier  in  the  thesis.  The  detection  program 
needs  to  be  made  more  robust  to  eliminate  the  event  time 
error.  A  scheme  Tor  ignoring  artifacts  should  be  included. 
The  location  program  wastes  time  looking  in  the  wrong 
direction,  although  the  bearing  accuracy  of  the  program  is 
very  goad.  The  algorithm  should  be  changed  to  quicklv  find 
the  right  bearing,  and  then  search  in  a  sector. 

The  type  of  each  event,  whether  it  was  a  pop  or  a 
whine,  was  not  recorded.  Collecting  this  information  and 
correlating  it  with  interarrival  time  and  range  still  needs 
to  be  done. 
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User's  Guide  -for  the  hdetect  Program 


Figure  A— 1:  Flow  Chart  of  the  hdetect  Program 


Source  Code  for  the  hdetect  Program 
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USER'S  GUIDE  FOR  THE  hdetect  PROGRAM 

The  purpose  of  the  hdetect  program  is  to  detect 
ambient  noise  transients  amidst  the  background  ambient 
noise  recorded  on  a  FRAM  data  tape.  This  is  done  bv 
comparing  the  short  average  of  data  points  to  the  long 
average  o-f  points  on  a  single  channel  in  order  to  flag  a 
possible  detection,  and  then  waiting  until  50V.  of  the 
channels  are  flagged  to  declare  an  actual  detection. 

The  input  for  the  hdetect  program  is  a  framread  output 
file  without  headers.  A  FRAM  data  tape  is  read  into  the 
file  by  the  command 

framread  -head  <  RETURN 

The  program  will  ask  for  the  input  device  (tape  drive 
designation) ,  the  output  file,  and  the  number  of  data 
segments  to  skip  and  to  read.  Each  segment  represents  ~.S 
seconds  of  data  on  24  channels.  The  framread  program  reads 
a  first  segment  which  contains  no  data  records,  so  you 
should  specify  skipping  one  more  segment  than  you  would 
normally  calculate.  For  example,  reading  the  entire  first 


half  of  a  20  minute  FRAM  IV  tape  would  require  the  response 
of 

1  160  < RETURN, 

to  the  Question  of  ''enter  #skip,  ^segments:". 

Once  this  input  file  has  been  created  the  hdetect 


l 


> 


program  can  be  used.  The  program  is  started  with  the 
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command 


hdetect 


< RETURN > 


The  program  will  ask  Tor  the  FRAM  tape  number,  the  Julian 
date  o-f  the  tape,  and  the  start  time  of  the  tape  in  hours, 
minutes  and  seconds.  The  program  will  then  ask  you  to 
select  the  channels  you  wish  to  use.  In  most  cases  the 
FRAM  data  tapes  did  not  have  ambient  noise  hydrophones  tied 
into  all  channels,  and  the  specific  channel  that  a 
hydrophone  was  recorded  on  changed  throughout  the 
experiment.  Which  channels  were  in  use  and  for  which 
hydrophones  can  be  found  in  the  experiment  logs.  The 
program  assumes  that  the  channel  number  is  equal  to  the 
hydrophone  number,  but  allows  you  to  change  this  by 
inputting  the  channel  number  and  the  proper  hydrophone 
number,  or  "O"  if  the  channel  is  not  in  use.  For  example, 
if  channel  3  was  not  used,  and  channel  7  was  used  for 
hydrophone  21,  the  input  would  be 


3,0 

7,21 

0,0 


< RETURN > 
< RETURN > 
< RETURN > 


The  "0,0"  ends  the  changes  to  the  channel  selection.  You 
must  now  hit  any  key  to  continue  the  program. 

You  will  be  asked  to  enter  the  input  device  (the  input 
framread  file) ,  the  number  of  skips  and  segments,  and  the 
name  of  the  output  file.  The  output  file  does  not  have  to 
exist  before  the  program  is  started.  It  will  be  created  by 
the  program.  The  number  of  skips  and  segments  are  those 


1 

ESI 


WTO 


that  you  Mould  calculate  using  3.8  seconds  per  segment. 

For  example,  to  proces  the  entire  first  half  of  a  20  minute 

FRAM  IV  tape  the  number  of  skips  and  segments  would  be 

0,160  < RETURN > 

This  is  all  of  the  input  required  by  the  user.  The 

program  proceeds  from  this  point  without  user  interaction. 

The  output  of  the  hdetect  program  is  a  file  containing 

a  list  of  detections  in  the  following  format: 

tapenumber  Juliandate  hour  minute  seconds 
0  eventnumber  eventtime 

channel  hydrophone  timedelay  amplitude 
channel  hydrophone  timedelay  amplitude 
channel  hydrophone  timedelay  amplitude 


channel  hydrophone  timedelay  amplitude 
0  eventnumber  eventtime 

channel  hydrophone  timedelay  amplitude 

channel  hydrophone  timedelay  amplitude 
0  eventnumber  eventtime 

-1 

The  "0"  at  the  start  of  a  line  indicates  a  new  event 
detection,  and  the  "-l”  at  the  start  of  a  line  indicates  an 
end  of  file.  Each  channel  that  was  flagged  for  a 
particular  event  is  listed  with  its  hydrophone  number, 
timedelay  from  the  earliest  channel  signal  arrival,  and  its 
peak  voltage  amplitude.  This  outfile  can  be  used  as  the 
input  file  for  the  location  programs  without  modification. 
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Figure  A-l  Flow  chart  o-f  the  hdetect  program 


/*  startdoc 
hdetect .  c 

program  to  read  whoi-segy  format  data  tapes  from  the  fram  IV  program, 
usage : 

hdetect 

program  is  interactive. 


by  Mary  Townsend-Manning 

enddoc 

*/ 


♦include  <atdio.h> 

♦include  <math.h> 

♦define  NCHAN  25 

♦define  RECLN  950  /*  number  of  samples  per  trace  (4  bytes  per  sample)*/ 

♦define  OBYTES  3800  /*  number  of  bytes  per  record  output  */ 

♦define  ZERO  0 

♦define  LONGFILTLN  64  /*  length  of  long  average  filter  */ 

♦define  FLN  64 
♦define  MAXE VENTS  4 
♦define  R£SET_DELAY  0.3 
♦define  RATIO  2.38 
♦define  3IGJ0ELAY  0.02 
♦define  THRESHOLD  0.5 
♦define  EVENT_DELAY  0.5 
♦define  END  -1 

double  vconv(x,y,n) 

register  floac  *x,  *v; 

register  int  n; 


( 

register  double  sum  -  0.; 


if (n  >  0) 

{ 

do  ( 

sum  +-  *x++  *  *y — ; 
)  while ( — n  >  0) ; 

} 

return (sum) ; 


main  ( ) 

i 
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float  l_ave [NCHAN] ,  sh_ave,  timeof flag [NCHAN] ; 
float  ampof flag (NCHAN] ,  firattime,  timedelay [NCHAN] ; 
float  chaneventtime [MAXEVENTS] (NCHAN] ; 

float  chaneventamp [MAXEVENTS ] [NCHAN],  eventime [MAXEVENTS ] ; 
int  flag (NCHAN] ; 
int  num_active_events; 
int  m,  n  ; 

int  event_f lag [MAXEVENTS 1 (NCHAN] ,  flag_aum,  number_of_events ; 
int  event_number (MAXEVENTS] ; 
int  nchan  “  0; 

int  chan [NCHAN] ,  k,  j,  i,  uid,  date,  hour,  min,  sec,  1; 
int  channel,  data,  toggle; 
char  answer; 

float  h[FLN],  longfilt [LONGFILTLN] ; 
int  nskip,nseg  ; 

char  oddobuf [NCHAN] [OBXTES] , evenobuf [NCHAN] [OBYTES] ; 
int  count, error; 
char  fname [80] , iname [80] ; 

FILE  *iptr,  *ptr,  *fp,  *fopen(); 
float  timeaeriea [NCHAN] [R£CLN+2*FLN] ; 
float  time  -  0.; 

num_active_eventa  -  0; 

for  (i"l; i<25; i++)  ( 
chan[i]-(i)  ; 
flag [i] -  0  ; 

for (m-l;m<MAXEVENTS;m++)  ( 
event_flag[m] [i]  -  0; 
chaneventtime [m] [i]  -  0; 
chaneventamp [m] [i]  -  0; 

1 

) 

/*  Program  initialization  from  keyboard  */ 


fprintf (stderr, "Program  Initialization\n" )  ; 
fprintf (stderr, "enter  FRAM  tape  »\n") 
f3canf  fatdin, "%d",  4tid)  ; 

fprintf (atderr, "enter  Julian  datein")  ; 
facanf (3tdin, "%d", sdate) ; 

fprintf (stderr, "enter  time  -  HR, MN, 3C\n" )  ; 
facanf (stdin, "%d, %d, %d", (hour, Smin, 4sec) ; 

fprintf  (stderr ,  "default  values  for  channels  and  phonesW); 
fprintf (atderr, "are  channel  #  «  phone  #.\n"); 

fprintf (stderr, "enter  channel,  phone  to  change. \n"); 
fprintf  (stderr,  "enter  '0'  for  phone,  to  eliminate  a  channel . \r.” )  ; 
fprintf (stderr, "enter  '0,0'  to  quit.Nn"); 
fscanf  (stdin,  "%d,  %d",  4  j,  4k)  ; 

while  (j  !-  0  SS  j  <  25)  ( 

chantj]  •  k; 


fscanf  (stdin,  "Id,  Id",  4  j,  4k)  ; 


fprintf (stderr, "FRAM  TAPE  Id  Julian  Date:  ld\n" , tid, date) ; 
fprintf (stderr, "  Time:  Id: %d: ld\n", hour, min, sec)  , 

for (i-1; i<13; i++) 

fprintf (stderr, "CH  Id  PH  Id  CH  Id  PH  ld\n", 

i, chan [i] , i+12, chan  I i+12) ) ; 

/*  Check  to  make  sure  inputs  are  correct  —  Change  if  necessary  */ 
fscanf (stdin, ”lc",4answer) ; 

fprintf  (stderr,  "Hit  any  key  and  RETURN,  when  ready."); 
fscanf (stdin, "Ic", {answer) ; 


for (i-1; i<NCHAN;i++)  ( 

if(chan(i]  !-  0)  nchan++; 


fpnntf  (stderr,  "enter  input  device:  ")  ; 
scanf ("Is", iname) ; 

if((iptr  -  f open ( iname,  "r"))  --  NULL) 

( 

fprintf (stderr, "can' t  open  ls\n", iname) 
exit ( 1 )  ; 


fprintf (stderr, "enter  #skip,  *segments:  \n")  ; 

fprintf (stderr, "values  of  0  and  320  will  read  entire  tape\n”); 
scanf ("Id, Id", 4nskip, 4nseg) ; 

/ *  load  bandwidth  filter  •/ 


if ( (fp-fopen ("PMfloat", "r") )  --  NULL)  ( 

printf  ("cannot  open  bandwidth  filter  file\n") 
exit (0); 

1 

for (i-0 ; i<64 ; i++) 

fscanf(fp,"lf",4h[i]) ; 
f close ( fp) ; 


/*  load  averaging  filter  */ 


for (i-0 ; i<LONGFILTLN; i++) 

longfilt'i]  -  1 . 0/ (float) LONGFILTLN; 


/*  Open  output  file  */ 


Source  code  for  the  hdetect 


program. 


•  u  » .» IT'i.'vi.'VYV' 


-124- 


fprintf (stderr, "enter  out-file:  ")  ; 
soanf ("%s", fname) ; 

if  (  (ptr  -  f open  (fname,  " w "))  --  NULL) 

( 

fpcintf (stderr, "can1 t  open  %s\n", fname) ; 
exit (1) ; 


fprintf (ptr, "%d  *d  td  *d  »d\n",  tid,  date,  hour,  min,  sec) 

time  -  time  +  3 . 8  *  (nsk.ip-1)  ; 

/*  enter  first  record  */ 

if  (nskip%2  —  1)  ( 

for ( j— 1 ; j<NCHAN; j++)  ( 

f read (ievenobuf ( j] (Q) ,  sizeof (float) ,  RECLN,  iptr); 

) 

toggle  -  1; 


else  { 

for( j-1; j<NCHAN; j++)  ( 

f read (ioddobuf ( j] [0] ,  sizeof (float)  ,  RECLN,  iptr); 

I 

toggle  -  0; 


/*  ENTERING  RECORD  READING  MODULE  */ 
time  *  time  +  0.504; 

/*■  fprintf  (stderr,  "using  buffer  size  %d  bytesXn" ,  sizeof  (buf!  ) 

f or ( i-1 ;  i  <  nseg;  i++) 
l 

/*  fprintf (stderr , "time  -  %f\n",  time); 

fprintf (stderr, "processing  record  *d\n",  nskip+i) ;  */ 

/*  read  next  record  into  appropriate  buffer  •/ 

if  vtoggle  —  1)  ( 

for (j-1; j<NCHAN; j++)  ( 

f read (4oddobuf ( j ] [0 ) ,  sizeof (float) ,  RECLN,  iptr); 

) 

toggle  -  0  ; 


else  ! 

for < j-1; j<NCHAN; j++)  ( 

f read ( ievenobuf ( j J (0 ] ,  sizeof ( float) ,  RECLN,  iptr) ; 
I 

toggle  -  1; 


/*  filter  and  square  data  */ 
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) 


if  (toggle  —  1)  ( 

for(j-l;  j<NCHAN;  j++>  { 

if(chan[j]  !-  0)  ( 

aq_f ilt (Soddobuf [ j] [Q] , Sevenobuf  ( j] [0] , h, 

Stimeseries [ j] [ 0 1 > ; 

l_ave(j]  -  vconv (Slongf ilt [0] , Stimeseries Cj] [2* (FLN-1) ) , 
LONGFILTLN) ; 


else  { 

for( j-1; j<NCHAN; j++)  ( 
if(chan[j]  !-  0)  ( 

sq_filt (Sevenobuf  C  j 3  [0] , Soddobuf  ( j  ]  [0)  ,  h, 

Stimeseries ( j]  (0] )  ; 

1  ave[j]  -  vconv  (Slongf  ilt  [01  ,  Stimeseries  [  j]  [2"  (FLN-1)  ]  , 
LONGFILTLN) ; 

) 

1 


/*  ENTERING  EVENT  DETECTION  MODULE  */ 

for  (k-0 ;  k<R£CLN;  1c  +-  4)  { 

for (1“1; 1<NCHAN; 1++)  ( 

if  (chan ( 1]  !-  0)  ( 

l_ave [I]  -  (63 . 0*l_ave [1]  +  timeseries [1] [2* (FLN-1) +k] ) /64 . 0; 
sh_ave  «  (timeseries  il)  [2*  (FLN-1)  +k]  + 
timeseries  [1]  [2*  (FLN-1)  +k-l]  +  timeseries  [1]  [2*  (FLN-1)  +k-2]  + 
timeseries  (1]  (2*(FLN-l)+k-3])/4.0; 

/*  reset  old  flags  '/ 

if  (flag (11  --  1  SS  (trme  -  timeofflagfl] )  >  RESET_DELAY) 
flag ( 1 J  -  0; 


/*  3et  flag  if  RATIO  of  signals  is  reached  */ 


if  ( <sb_ave/l_ave[l))>«RATIO)  ( 
if  (flag(l)  —  1)  ( 

if  (sh_ave  >  ampof f lag [ 1J )  ( 

timeofflagtl]  -  time; 
ampofflagfl)  -  sh_ave; 

( 

) 

else  ( 

if  (num_act ive_e vents  — —  0)  ( 

flag ! 1]  -  1; 
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timeoff lag [1]  “  time; 
ampofflagtl]  _  sh_ave; 

) 

else  { 

for  (m-l;m<  (num_active_e vents  +  1) ;m++)  ( 

if  (event_f  lag [m]  [1]  "  1)  ( 

if  {(time  -  chaneventtine [m) [1] ) 

<-  SIG_DELAY)  { 

if  (sh_ave  >  chaneventamp [m] [1])  t 
chaneventtime tm] [1]  “  time; 
chaneventamp [m] (1)  -  sh_ave; 

) 

1 

else  { 

flag  [1]  -  1; 

timeof flag [ 1]  ■  time; 

ampofflagtl]  ”  sh_ave; 


else  { 


event_f  lag [m]  [1]  _  1; 
chaneventtime [m] [1]  ■  time; 
chaneventamp [m] [1]  “  sh_ave; 
m  -  num  active  events; 


/*  end  of  set  flag  module  */ 
/*  start  new  event  module  */ 
f lag_sum  “  0 ; 


for  (1-1;1<NCHAN;1++)  ( 

if  (flag [1]  --  1)  flag_sum++; 

if  (( (float) flag_sum/ (float) nchan)  >-  THRESHOLD)  ( 
num_active_events++; 

eventime lnum_active_events]  “  time; 
number_of_events++; 

event_number (num_active_events ]  -  number_of_events ; 

for  U-1;1<NCHAN,-1++)  { 

event_flaglnum_active_events] 11]  ■  flagll]; 
chaneventtime [num_active_events]  11]  “  timeof flag  11] 
chaneventamp [num_active_events] [1]  ”  ampofflagtl]; 
flagll]  -  0; 
timeof flag [ 1]  ”  0; 
ampofflag(l)  ”  0; 
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/*  end  of  new  event  module  */ 

/*  start  of  deactivate  old  event  module  */ 

if  (num_active_evencs  >  0  66  {time  -  eventimefl])  >  EVENT_DELAY)  { 
fprintf  (ptr,"ld  Id  If \n", ZERO, 

e vent_n umber (1] , event ime [1) ) ; 

/*  fprintf  (stderr,”ld  If \n" , event_n umber ( 1] , event ime [1] ) ;  */ 


find  time  delays  by  finding  earliest  channel  event 
time,  and  subtracting  that  from  the  other  channel 
times 


firsttime  -  10000.0; 

for  (1-1;  1<NCHAN;  1++)  { 

if(chan(l]  !-  0  66  event_f  lag  1 1]  [  1]  !-  0  66 

chaneventtime ( 1 )[ 1]  <  firsttime) 
firsttime  -  chaneventtime  1 1]  ( 1]  ; 


for  (1-1;  KNCHAN;  1  +  *)  ( 

if(chan(l]  1-  0  66  event_f lag  1 1 ]  ( 1)  !-  0)  I 

timedeiaytl]  -  (  (chaneventtime [1]  [ 1] )  -  firsttime), 
fprintf (ptr, "Id  Id  If  If \n” , 1 , chan ( 1 ) , 
t imedelay [ 1] ,  chane vent amp  11] 1 1 j ) ; 


print  to  file  to  indicate  end  of  event  */ 

for  ( 1-1 ; l<num_active_events ; 1++)  1 

for  (m-1;  m<NCHAN;m++)  { 

event_f lag [ 1] (m)  -  event_f lag [ 1+1 J [m] ; 
chaneventt ime { 1 J [m]  -  chaneventoime [1+1 ] [m] , 
chaneventamp  [1]  [m]  -  cr.aneventamp  1 1  +  1]  [m]  ; 

I 

eventime [1]  -  eventime [1+1] ; 

even t_n umber [1]  -  event_number [1+1] ; 

) 

num_active_e vents — ; 

) 

t ime  +—  0.004*4.0; 
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/*  fprintf (stderr,  "processing  record  %d\n",  nskip+nseg) ;  */ 
if (toggle  ~  1)  ( 

for(j-l; j<NCHAN; j++)  ( 

if(chan(j]  !-  0)  ( 

for (k-0;k<OBYTES;  k++) 
oddobuf [jj (k)  -  0; 

sq_f ilt (tevenobuf [  j  ] [0 ] , t oddobuf [  j  ] [0] , h, ttimeseries [  j  ] [0 ] ) ; 

l_ave(j]  -  vconv (tlongf ilt [0] ,  ttimeseries [ j] [2* (FLN-1) ) , 
LONGFILTLN) ; 


for  (  j-1;  j<NCHAN;  jm+)  ( 
if(chan[j]  ! -  0)  ( 

for (k— 0 ; k<OBYTES ; k++) 
evenobuf ( j] [k]  -  0: 

sd_f  ilt (t oddobuf [ j]  [0] , tevenobuf t  j  ]  [0 ] , h, ttimeseries ( j  J  [0  3 ) ; 

l_ave(j]  -  vconv (tlongf ilt [ 0 ] , ttimeseries  \  j ]  [2* (FLN-I) 3 , 
LONGFILTLN) ; 


for (k-0; k< (RECLN-2* (FLN-1) > ;k+-4>  ( 

for  ( 1-1;  KNCHAN;  !■>■*)  ( 

if  (chan[i;  0)  ( 

i_ave(I3  ~  i 63 . 0 * l_ave ! 1 ]  *  time series ;  1  ]  ; 2 *  r FLN- 1 ) -< ' v 
sh_ave  -  ttimeseries ( 1 ) (2* (FLN-1 ) *k ]  * 
timeseries  [1]  (2*  (FLN-1) +k-l]  *  timesenes  ( 1]  (2*  (FLN-1)  *k-2  ;  - 
timeseries [1) 32* (FLN-1) +k-3) ) /4.0; 

/*  reset  old  flags  */ 


(fiagii;  --  1  it  (time  - 
flag  Ilj  -  0; 


timeof flag  ]  1 ]  >  ?E2ZT  ce; 


flag  if  RATIO  of  signals  is  reacned  * 

if  ( (sh_ave/l_ave [ 1] ) >-RATlO)  I 
if  (flagdj  --  1)  { 

if  (sh_ave  >  ampof f lag ; 1 3 )  ( 

timeof f lag i 1 i  -  time; 
ampofflag;!;  -  sh_ave; 


if  (num_active_events  0) 
flag;!]  -  1; 
timeof flag ( 1'  -  time; 
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ampofflag[ll  -  sh_ave; 


) 

else  ( 

for  (m-1 ;m< (num_active_events  +  1) ;m++)  i 
ii  (event_f lag (m) [ 1]  —  1)  ( 

iS  ((time  -  chaneventtime [m] [ 1) ) 

<-  SIG_DELA¥>  ! 

if  (sh_ava  >  chaneventamp (ml [1] ) • 
chaneventtime (m) [1 ]  -  time; 
chaneventamp (m] [1]  -  sh_ave; 

I 

I 

else  ! 

flag [1]  -  1; 
t imeof f lag  [  1]  -  time; 
ampofflag[l]  -  sh_ave; 

) 


else 


event_f lag [m] ( 11  -  1; 

chaneventtime fm] (1)  -  time; 
chaneventamp [m] [1]  -  sh_ave; 
m  -  nun  active  events; 


1 


I 


end  of  set  flag  module  */ 
3tart  new  event  module  */ 


f lag_sum 


0; 


for  (1-1;  1<NCHAN;  1+M  ( 

if  (fiag(l)  1)  fiag_sum++; 

1 

if  (( (float) flag_sum/ (float) nchan)  >-  THRESHOLD)  ( 
n  um_a  c  t  i  ve  _e  ve  n  t  s  +  * ; 
eventime (num_active_eventsl  -  time; 
number_o  f _e vent s ++ ; 

event_number [num_active_events]  -  number_of_e vents ; 

for  (1-1;  KNCHAN;  1++)  ( 

event_flag(num_active_events) (11  “  flag(l); 
chaneventtime  [num_active__events]  [1)  -  timeofflagl 
chaneventamp (num  active_events) ( 1)  -  ampofflag(l) 
Slag [1]  -  0; 
timeof f lag ( 1]  -  0; 
ampofflag[l]  -  0; 
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/*  end  of  new  event  module  */ 

/*  start  of  deactivate  old  event  module  */ 

if  (num_active_e vents  >  0  44  (time  -  eventime (1])  >  EVENT_DE1AY)  ( 
fprintf  (ptr,"%d  td  %f\n",  ZERO, 

event_number  £  1] , eventime  [  1] )  ; 

/*  fprintf  (stderr,"%d  %f\n” , event_number [1] ,  event  ime [ 1 ]) ;  */ 


find  time  delays  by  finding  earliest  channel  event 
time,  and  subtracting  that  from  the  other  channel 
times 


firsttime  -  10000.0; 

for(l-l;l<NCHAN;l++)  ( 

if  (chan  (1]  0  ti  event_f  lag  [  1)  [  1] 

chaneventtime t 1] [1]  <  firsttime) 
firsttime  “  chaneventtime (1)[1J; 


for  ( i”l ;  1  <NCHAN;  1  **)  ( 

if (chan [1]  !-  0  4S  event_flag ( 1] ; 1]  0)  ( 

timedelay(l)  «  (  (chaneventtime [ 1 ]  ti] )  -  firsttime) 
fprintf (ptr, "%d  %d  %f  %f \n" , 1, chan  ( 1] , 
timedelay { 1 ) ,  chaneventamp [ 1 J [ 1 ] ) ; 


print  to  file  to  indicate  end  of  event  •/ 

for  ( 1-1 ; l<num_active_events ; 1+*)  i 
for  (m—1;  m<NCHAN;mf-»)  ( 

event_f lag ( 1] [m]  -  e vent_f lag ! 1*1 ] [m] ; 
chaneventtime ( 1 ] (m]  -  chaneventtime C i*l ] (m! ; 
chaneventamp ( 1] (m)  -  chaneventamp [ 1+1 ] (m] ; 

! 

eventime(l)  -  eventime [ 1*1 ] ; 
event_numbec ( 1]  «  event  number [ 1*1] ; 

I 

num  active  events--; 


time  +"  0.004*4.0; 


/*  print  out  all  events  */ 
if (num  active  events  >  0)  I 
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for (k-l;k< (num_active_events  +  l);k++)  ( 

fprintf  (ptr, "td  id  %f\n",  ZERO, 

event_number  [k]  ,  eventime  [k] )  ; 

/*  fprintf  (3tderr,"%d  %f\n" , event_number [k] , eventime [k] ) ;  */ 


find  time  delays  by  finding  earliest  channel  event 
time,  and  subtracting  that  from  the  other  channel 
times 


firsttime  «  10000.0; 

for  (1-1;  KNCHAN;  1++)  1 

if(chan[l]  !-  0  is  event_f  lag  (k]  f  1  ]  !-  0  it 
chaneventtime t k]t 1]  <  firsttime) 
firsttime  »  chaneventtime [ k ] [11; 


for  (1-1;  KNCHAN;  1++)  ( 

if(chan(l]  !-  0  Si  event_f  lag  [k)[  1]  !-  0)  ( 

timedelay ( 1]  »  {  (chaneventtime [k] 1 13 )  -  firsttime) 
fprintf  (ptr,  "Id  id  %f  If \n", 1, chan [1] , 
timedelay [ 1) ,  chaneventamp (k] ( 1] ) ; 


/*  final  summary  to  screen  */ 


fprintf (ptr, "%d" 
fclose(ptr)  ; 
fclose(iptr)  ; 
exit (0)  ; 


sq_f ilt (first , second, filter, output ! 

/*  filters  and  squares  two  data  arrays  */ 


float  first [RECLN] ,  second [ RECLN] ,  filter[FLNJ 
float  output ( AECLN+2*FLN) ; 


Source  code  -far  the  hdetect  proaram. 


L>  Tji  -JO\*  -JT’JfJ-j. -Jl-J  -> 


.  'T.  lT.  »■. 


.  •  ■  'T. 
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/*  put  zeros  in  first  (FLN  -1)  places  of  output  */ 
int  j,  i; 

float  transition [2 *FLN] ; 
float  sum; 

for (i-0 : i<FLN-l ; i++) 

output [i]  “0.0  ; 

/*  put  in  first  data  */ 

for(i-FLN-l;i<RECLN;i++)  ( 

sum  “  vconv  (sf  ilter  [01 , 4f  irst  [i] ,  FLN)  ; 
output [i]  ”  sqrt (sum  *  sum) ; 

1 

/*  put  in  transition  from  first  to  second  data  */ 

for ( j~0; j<FLN-l; j++)  1 

transition  ( j]  •  f irst (RECLN- (FLN-1) + j 1 ; 
transition ( j+FLN-1]  -  second  [j].* 


for (i-0;i<FU»-l;i++) ( 

sum  “  vconv (sfilter [0 ] , Stransition [i+FLN-1] , FLN) 
output [i+RECLN]  -  sqrt (sum  *  sum) ; 

I 

/*  put  in  second  data  */ 

for (i-FLN-l;i<2*FLN;i++)  [ 

sum  -  vconv  (4f ilter [0] , 4second[ i] , FLN) ; 
output [RECLN+i]  “  sqrt (sum  *  sum); 


Source  code  for  the  hdetect  program. 


USER'S  GUIDE  FOR  THE  location ,  far  locate , 


and  fine locate  PROGRAMS 

The  purpose  of  the  location  programs  is  to  -find  the 
spatial  location  of  an  event  -from  the  time  delays  between 
signal  arrival  at  different  hydrophones.  The  program 
assumes  a  test  location  and  computes  the  slant  range  to  the 
individual  hydrophones.  The  slant  ranges  are  platted 
against  the  experimental  time  delays  and  a  least  squares 
fit  is  done.  The  test  location  with  the  best  least  squares 
fit  is  considered  the  location  of  the  event. 


The  input  to  the  location  programs  is  the  output  file 
of  the  detection  program.  Manual  time  delays  may  be 
substituted  for  the  program  generated  time  delays  in  this 
file,  but  this  editting  must  be  done  before  the  location 
program  is  invoked. 

This  location  program  is  very  user  interactive.  The 
user  starts  the  location  program  with  the  command 
location  < RETURN > 

The  program  asks  for  the  input  and  output  file  names.  It 
then  reads  the  input  file  and  asks  whether  the  user  would 
like  to  locate  the  first  event.  This  allows  the  user  to 
skip  down  to  the  event  of  interest.  The  program  then 
allows  the  user  to  adjust  which  hydrophone  time  delays  wi 1 l 
be  used  in  the  location  process.  This  is  very  hand,  for 
removing  questionable  time  delays,  in  order  to  get  a  better 
location  rdution.  U:  ►h  the  hvdroDhone  channels  rhnspn. 


the  program  proceeds  with  the  actual  location  algorithm. 


In  the  fznelcoate  and  far  locate  programs  the  user  is  asked 
to  specify  which  quadrant  or  direction  is  to  be  searched. 

The  location  program  tries  test  locations  in  a  large 
grid,  and  when  the  "best"  location  is  -found,  then  searches 
a  smaller  grid  around  this  "best"  location.  The  location 
and  far  locate  programs  have  4  levels  o-f  grids  and  the 
fmelocate  program  has  3  levels.  Intermediate  answers  are 
displayed  tor  each  level. 

The  intermediate  and  -final  answers  display  the  and  y 
coordinates  of  the  best  location,  the  standard  deviation  of 
the  least  squares  fit  (sigma),  the  group  speed  (which 
should  be  around  1440  m/si .  and  the  /  interceDt  of  the  time 
delay  /  slant  range  plot.  After  the  final  answer  the  user 
is  asked  whether  or  not  he  would  like  to  remove  outlying 
points.  If  this  action  is  selected  the  program  removes 
hvdrophones  with  a  deviation  from  the  least  squares  fit  of 
more  than  3  times  the  standard  deviation,  and  the  group 
speed  and  the  standard  deviation  are  recalculated  and 
di spl ayed . 

The  user  is  then  asked  if  he  would  like  to  locate  the 
event  with  different  hydrophones,  and  if  so  returns  the 

user  to  the  start  of  the  channel  selection  process.  The 
user  m  a  .  rep  eat  tnis  location  scneme  as  mane  times  as 
necessary  for  a  particular  event.  Once  the  user  is 
■a*-’  a  f  j  p».-|  i-i  i  tf-.e  1  oc  1 1  on  -n  =  wer  .  and  deal  *  ne«.  to  l  a- -at  a 
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the  event  with  different  phones,  the  program  calculates  the 
event  strength  based  on  the  amplitudes  in  the  input  file, 
the  location  of  the  event  and  spherical  spreading  lasses. 
The  location  parameters  and  the  strength  are  then  written 
to  an  outfile. 

At  this  point  the  user  is  given  the  option  to  exit  the 
program  or  locate  the  next  event  in  the  input  file.  The 
location  process  continues  until  the  user  exits  or  until 


the  last  event  in  a  file  is  located. 
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/*  location. c 

program  to  locate  the  spatial  position  of  an  event  based  on 
time  delays  taken  from  "detection" . 

Source  strength  is  also  computed. 

*/ 

♦include  <stdio.h> 

♦include  <math.h> 

♦define  PHONES  31 
♦define  LEVEL  4 
♦define  FINENESS  20 
♦define  DEPTH  91.0 

♦define  SENSITIVITY  0.0000000112202 
main  ( ) 


float  amptPHONESl,  phonex (PHONES; ,  phoney [PHONES ) ; 

float  timedelay [PHONES] ,  r [PHONES],  bestrange [PHONES ) ; 

float  sumtime,  sumr,  sumrsq,  sumtimer,  slope,  yintr; 

float  sigma,  bestsigma,  bestslope,  bestyintr; 

float  bestamp,  N,  source [PHONES ] ,  sumsource; 

float  xgs,  yga,  xcntr,  ycr.tr,  a,  b,  bestx,  besty; 

float  level,  gridsize,  xfineness,  yfineness,  time,  gpspeed; 

int  i,  j,  tape,  date,  hour,  min,  sec,  event,  flag; 
int  phone,  num,  1,  n,  m,  bestflag; 
int  phonef lag [PHONES )  ; 

int  eventselect,  rerun,  change,  answer,  bye,  quadrant; 
char  iname[80],  oname[80J; 

FILE  *ptr,  open O ,  "locptr,  'optr; 

/*  PROGRAM  INITIALIZATION  */ 
bye  »  0 ; 

fprintf (stderr,  "input  file  -  \n")  ; 
scanf("%s",  iname); 
fprintf (stderr, "output  file  ■  \n")  ; 
scanf ("\s",  oname) ; 


/*  open  files  */ 

if ( (ptr  -  fopen (iname,  "r”) )  ~  NULL)  ( 
fprintf (stderr, "can' t  open  %s\n”,  iname); 
exit  (1)  ; 

! 


Source  code  for  the  far  locate  program 


V  *  *  '  rfv  .  -  r.  «.  rt  v-_  r.  «r,  r.  *,  ,  .  f,  r» 


i 

A 
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if ( (optr  -  fopen(oname,  "w"))  --  NULL)  t 
fprintf  (atderr,  "can' t  open  %a\n",  oname) ; 
exit  (1)  ; 

) 

if((locptr  -  fopen ("array_loc",  ”r"))  --  NULL)  1 
fprintf (atderr, "can' t  open  array_loc  file\n") ; 
exit ( 1 )  ; 

) 


/*  read  hydrophone  locations  into  array  */ 
for  (i-1;  i<PHONES;  i++)  { 

f3canf (locptr, ”%d  %f  %f",Sphone,  Sphonex[i],  Sphoney  [i] ) ; 

) 


/*  read  input 
fscanf (ptr, 
/*  read  event 
f acanf (ptr, 


file  header  */ 
"%d  %d  %d  %d  %d", 
header  */ 

"%d  ",  sf lag) ; 


Stape,  Sdate,  Shour,  imin,  Ssec)  ; 


while (bye  !-  -1)  1 


£m 


'i 


rm 

J 


i 


eventselect  -  0; 

while (eventaelect  !-  1)  { 

if (flag  <  0) 
exit (0) ; 

for  (i-1;  i<PHONES ;  i++>  ( 

phoneflagti]  -  0; 

) 

fscanf (ptr,  "%d  %f",  ievent,  itime); 

for  (i-1;  KPHONES;  i++)  ( 

fscanf (ptr, "%d", Sflag) ; 


V 
>  * 
> 

V 
/ 

V 

V 


if (flag  >  0)  (  v 

fscanf  (ptr,  "»d",  Sj); 
phoneflaglj]  -  1; 

fscanf  (ptr,  "*f  %f",  stimedelay [ j] ,  Samp(jJ); 

1  S 

else  ( 

i  -  PHONES; 

I 


fprintf  (stderr,  "event  -  %d,  time  -  %f\n",  event,  time); 

/. 

r. 

/ 


Source  code  for  the  far  locate  program 
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fprincf (stderr, "Do  you  wish  to  locate  this  event?  (1  -  yes)\n"), 
scanf ("%d",  4eventselect) ; 


/*  channel  selection  */ 


rerun  m  1; 

while  (rerun  1)  ( 

change  -  1; 

while  (change  “**  1 )  1 


fprintf (stderr, "  phone 
for  (i“l;  i<PHONES;  i++)  ( 

if (phoneflag (i]  1-0)  ( 


fprintf (stderr, 


%f  \n" , i, timedelay [ ij  )  ; 


fprintf (stderr, "Do  you  wish  to  change  status?  (1  •  yes)\n"); 
scanf ("%d",  ichange) ; 
if  (change  --  1)  ( 


fprintf (stderr, "change  status  by  typing  phone»\n"); 
fprintf (stderr, "type  -1  to  quit\n") ; 
scanf ("%d", 4  j )  ; 


while  (j  !“  -1)  ( 

if  (phoneflag!  j]  !-  0)  f 
phoneflag i;]  “  0; 


else  ( 

phoneflagtj]  -  1; 


3canf ( "%d",  4  j  ]  ; 


/*  locate  event  */ 


num  ■  0  ; 
sumtime  -  0.0; 


for  (i-1;  i<PHONES;  i++)  ( 

if  (phoneflag  i  i]  !-  0)  ( 


sumtime  +-  timedelay li) ; 


fprintf (stderr, "select  quadrant  to  search  (1-NE,  2-NW,  3-SW,  4-SE\n") , 
fprintf (stderr, "5-N,  6-S,  7-E,  8«W)\n"); 
scanf ("%d" , 4quadrant) ; 
if  (quadrant  --  1)  ( 


Source  code  for  the  far  locate  proqram 


»»v  ^ 


' 

' 

r 

fc 

> 


l 


| 


i 

i 


n  ir»  k-«  vvv.’v 


.v.vv.wv  v  *.* v.*  K*  WJ *T(.^LW.7W.*V.*V 
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xcnt  r  -  10000.0; 
ycntr  ”  10000.0; 

) 

else  if  (quadrant  — 
xcntr  “  -10000.0; 
ycntr  »  10000.0; 

! 

else  if  (quadrant  — 
xcntr  "  -10000.0; 

ycntr  ■  -10000.0; 

1 

else  if  (quadrant  - 
xcntr  -  10000 . 0 ; 
ycntr  —  -10000.0; 


else  if  (quadrant  — 
xcntr  «  0.0; 
ycntr  —  10000.0; 

else  if  (quadrant  —  s> 
xcntr  -  0.0; 
ycntr  -  -10000.0; 

I 

else  if  (quadrant  —  ') 
xcntr  -  10000.0; 
ycntr  -  0.0; 

else  if  (quadrant  —  8) 
xcntr  -  -10000.0; 
ycntr  “  0.0; 

! 


else  ( 

xcntr  ”  0.0; 
ycntr  -  0.0; 

) 

bestflag  “  0; 

for (1-0;  1< LEVEL;  1++)  l 
level  —  LEVEL— 1—1. 
gridsize  —  pow (10 . 0, level)  . 

for (m-0 ;  m<FINENESS;  m++)  ( 

yfineness  -  m  -  (FINENESS/2) ; 

ygs  “  ycntr  +  yfineness  *  gridsize; 

for (n-0;  n<FINENESS;  n++)  ( 

xf ineness  -  n  -  (FINENESS/2), 
xgs  -  xcntr  +  xfineness  *  gridsize; 
sumr  -  0.0; 
suinrsq  -  0.0; 
sumtimer  -  0.0; 
for  (i-1;  i<P HONES;  i++)  < 

if (phonef lag [i]  !”  0)  ( 

a  -  xg3-phonex[i] ; 


Source  code  -for 


the  far  locate  program 
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b  “  ygs-phoney  fi]  <" 

r[i]  “  sqrt (pow  (a, 2 . 0)  +  pow(b,2.0) 


pow (DEPTH, 2.0) ) ; 


sumr  +-  r[i]; 

sumrsq  +—  pow  (r  [i] , 2 . 0)  ; 

sumtimer  +-  <t[i]  *  t ' medelay [i] ) ; 


N  -  num;  ,  .  , 

slope  -  ( (N  *  sumtimer)  -  (sumr  *  sumtime))/ 
( (N  *  sumrsq)  -  pow (3umr, 2 . 0) ) ; 
yintr  -  (sumtime  -  (slope  *  sumr) ) /N; 


s igma  -  0.0; 

for  (i-1;  KPHONES;  i++>  ( 
if (phonef lag (i]  !-  Q)  1 

sigma  +-  pow ( (timedelay til -yintr- (slope 

) 

I 

sigma  -  sqrt (sigma/N) ; 


r [ i 1 ) ) ,2.0) 


if  (bestflag  --  0)  l 
bestx  -  xgs; 
besty  -  ygs; 
bestsigma  -  sigma; 
bestslope  -  slope; 
bestyintr  -  yintr; 
for  (i-1;  KPHONES;  i++)  ( 

if  (phonef lag  ’.”0)  ( 

bestrange(i)  “  r[i]; 

I 

) 

bestflag  -  1; 


else  ( 

if  (sigma  <  bestsigma)  { 
bestx  -  xgs; 
besty  “  ygs; 
bestsigma  -  sigma; 
bestslope  -  slope; 
bestyintr  *  yintr; 
for  (i-1;  KPHONES;  i++>  I 
if  (phoneflag  !-  0)  ( 

bestrange[i]  -  r(il; 

I 

) 


) 


Source  code  for  the  farlocate  program 


wu  n  <rv  wv,  1 r.y.  w  «r»  v,  *W.  <r.  I'.  *. 


gpspeed  -  1 . O/bestslope ; 


fprintf ( Stderr, "bestx  -  \t ,  besty  -  If,  sigma  -  %f\n", 
bestx,  besty,  bestsigma) ; 
fprintf ( stderr, "group  velocity  «  %f\n",  gpspeed); 
fprintf (stdecr, "y  intercept  -  %£\n",  bestyintr); 
xcntr  -  bestx; 
ycntr  -  besty; 


f printf ( scderr ,  "Oo  you  wish  to  remove  outlying  points?  (1-yes) \n”); 
scanf ("%d”, Sanswer) ; 
if  (answer  --  1)  ( 

num  -  0 ; 
sumr  -0  . ; 
sumt  ime  -  0  .  ; 
sumrsq  -  0 . ; 
sumtimer  -  0 .  ; 

for  (i-1;  KPHONES;  i++)  ( 

if  (phonef lag [i]  !-  0)  ( 

if  (sqrt (pow ( (timedelay (i) -bestyintr- (bestslope'bestrange [ i] ) ) 
2.0))  <  2 . 5*bestsigma)  l 
num+*; 

sumtime  +-  timedelay  [i!  ; 

sumr  +-  best  range ti] ; 

sumrsq  best range  I i ) ‘best  range t i] ; 

sumcimer  ■*—  best  range  [  i  ]  ‘timedelay  [ij; 


else  ( 

fprintf (stderr, "outlying  phone  *  %d\n",i); 


N  -  num; 

bestslope  -  ( (N'sumtimer) - (sumr* sumt ime) ) / 

( (N*sumrsq) - (sumr'sumr) ) ; 
bestyintr  -  (sumtime- (bestslope'sumr)  ) /N; 
gpspeed  «  1 . O/bestslope ; 

fprintf (stderr, "bestx  -  %f,  besty  -  *f,  sigma  -  %£\n" 
bestx,  besty,  bestsigma); 
fprintf (stderr, "group  velocity  -  %f\n",  gpspeed); 
fprintf (stderr, "y  intercept  -  Sf\n",  bestyintr); 


fprintf (stderr, 

"Do  you  wish  to  relocate  with  different  phones?  (1-ves) \n") 
scanf ("%d", s rerun) ; 


Source  code  for  the  far  locate  program 


vvv.v,. 


/*  finding  source  amplitude  */ 

for (i-1;  i<P HONES ;  i++)  ( 

if  (phonef lag [i]  !»  0)  ( 

scurce[i]  -  (amp  [ i) /SENSITIVITY)  *  bestrange [i] : 
sumaource  +”  source [i]; 

} 

) 

bestamp  -  sumsource/N; 

fprintf (optr, "%d  %d  \t  %f  %f  %f\n”,  event,  date,  time, 
beatx,  besty,  bestamp); 

fprintf (optr, "%f  %f  %f\n",  bescsigma,  gpspeed,  bestyintr); 

f  printf  ( stderr,  ”Oo  you  wish  to  exit’  (-1  -  exit)\n''); 
scanf("%d",  sbye) ; 

) 

f  close  (ptr)  ; 
f close (optr) ; 
f close ( locptr) ; 

exit (0)  ; 

) 


source  cade  far  the  Tar  ]  if.  a  t  e  program 
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APPENDIX  D 


Table  D— 1:  Angles,  Ranges  and  Times  for  Refractive 

Propagation  Paths 


'able  D-2:  Spreading  Loss  Function,  G(r) 
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Table  D-2  Spreading  Loss  Function,  Sir) 

Theta  zero  Theta  1  z  R1  R2  6(R1)  6IR2)  10  log  10  log 

G(R1)  6!R2) 

Theta  zero  is  the  surtace  launch  angle. 

Theta  1  is  the  angle  oF  the  retractive  path  at  the  hydrophone  depth  (93«> 
z  is  the  laxuua  depth  at  the  retractive  path, 

R1  is  the  horizontal  range  to  the  hydrophone  intersected  on  the  downward  swing 
R2  is  the  horizontal  range  to  the  hydrophone  intersected  on  the  upward  swing 
S(R)  is  the  spreading  loss  Function. 
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